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Resumen 


El grafeno es un cristal bidimensional consistente en una red hexagonal de átomos 
de carbono. Su estructura de bandas a baja energía y en situación de medio llenado 
viene descrita por la ecuación de Dirac sin masa en (2+1) dimensiones. Este hecho 
hace que gran parte de las propiedades electrónicas y de transporte que presenta este 
sistema sean diferentes a aquellas que se encuentran en los semiconductores de baja 
dimensionalidad usuales. 

En la realidad, las muestras de grafeno se encuentran curvadas. Las principales 
fuentes de curvatura son las fluctuaciones térmicas, los defectos topológicos y las 
interacciones con el sustrato. El hecho de encontrarse curvado hace que el grafeno 
presente modificaciones en sus propiedades físicas con respecto a la situación ideal. 

En esta tesis proponemos un modelo de fermiones de Dirac en espacios curvos para 
modelizar estas desviaciones devidas a la curvatura en las propiedades electrónicas y 
de transporte del grafeno. Hemos calculado perturbativamente la densidad de esta- 
dos local en presencia de defectos topológicos y en presencia de corrugaciones suaves. 
También hemos calculado la conductividad DC a nivel semiclásico en el grafeno cur- 
vado por defectos topológicos, encontrando un comportamiento de tipo difusivo, en 
contraste con la conductividad hallada para otros tipos de desorden estudiados en la 
literatura. 

Por último, proponemos una forma alternativa a las expuestas hasta ahora en la 
literatura para la obtención de haces electrónicos polarizados en el índice de valle, 
por medio del uso de la distorsión trigonal en las bandas del grafeno. Este efecto de 
separación de haces puede observarse en uniones n—p—n” bajo unas condiciones ac- 
cesibles experimentalmente. Por último, describiremos este mismo efecto en cristales 
fotónicos. 
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Capítulo 1 


Estructura de la tesis. 


El capítulo segundo consiste en una introducción al material llamado grafeno. 
Pondremos este material en el contexto de todas las formas alotrópicas del carbono. 
Revisaremos la manera en la que se modeliza la estructura electrónica de baja energía 
de este material a medio llenado, encontrando la muy interesante propiedad de que 
las excitaciones electrónicas a baja energía vienen descritas por la ecuación de Dirac 
sin masa en (2+ 1) dimensiones. La dimensionalidad del sistema y la forma cónica de 
la relación de dispersión E(k) se traducen en una densidad de estados que depende 
linealmente con la energía, siendo nula en los llamados puntos de Dirac (puntos en 
los cuales los vértices cónicos de las bandas de conducción y de valencia se tocan). 
La densidad de estados nula está en el corazón de la mayor parte de las propiedades 
anómalas que se observan en el grafeno. De entre todas esas propiedades anómalas 
desde el punto de vista de la teoría de semiconductores, la más interesante para 
nosotros será la respuesta que presenta el grafeno al desorden. 

En el capítulo tres nos introducimos en el estudio principal de esta tesis: el es- 
tudio de la estructura electrónica del grafeno en presencia de desorden topológico. 
Como sabemos del capítulo anterior, el grafeno consiste en una lámina bidimen- 
sional de átomos de carbono dispuestos en estructura hexagonal. Un defecto to- 
pológico consiste en la introducción de un polígono de lado distinto a seis en dicha 
estructura hexagonal. Los defectos más comunes son los pentágonos y heptágonos 
(también conocidos como disclinaciones), o defectos compuestos, como dislocaciones 
(un par heptágono-pentágono) o los célebres defectos Stone-Wales (doble pareja de 
pentágonos y heptágonos). 

Comenzaremos el capítulo respondiendo a la pregunta de si los cristales bidimen- 
sionales pueden existir. La respuesta es sí, gracias a la interacción entre los modos 
fonónicos flexurales y los modos fonónicos contenidos en el plano, pero pagando como 
precio el uso de la tercera dimensión, es decir, efectivamente, los cristales bidimensio- 
nales suspendidos son estables, pero presentarán fluctuaciones estructurales en forma 
de corrugaciones. Aparte, gracias a este doblado en la tercera dimensión, es posible 
argumentar que los defectos topológicos pueden existir, y por tanto, son una fuente 
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natural de curvatura para la lámina real de grafeno, junto con las propias fluctuaciones 
citadas antes. 

Tras la discusión anterior, nos metemos de lleno en el estudio de las consecuencias 
que tienen estos defectos sobre la estructura electrónica. El primer modelo que vamos 
a considerar consiste en la consideración de las condiciones de contorno sobre las 
funciones de onda generadas por estos defectos. Veremos que estas condiciones de 
contorno se pueden transformar en tres potenciales gauge ficticios. El primero está 
relacionado con la frustración de las subredes que se produce al introducir un polígono 
de lado impar El flujo “magnético” asociado a este campo es 7/2Rp, donde Rp es 
la matriz 72 % 02, la estructura peculiar de esta matriz tiene como consecuencia que 
los puntos de Dirac se mezclen. El segundo campo está relacionado con la falta o el 
exceso angular en el material producido al introducir un defecto. Los electrones en el 
grafeno son sensibles a dicho defecto o exceso de material, pudiendo entenderse este 
campo como una “modificación” en su momento angular, al estar asociado el flujo a 
la matriz 73 2 03 de rotaciones. El tercer campo gauge aparece cuando se consideran 
múltiples defectos, y su origen es la no conmutatividad de los operadores de transporte 
paralelo de la función de onda a lo largo de rotaciones en torno a cada defecto y de 
las traslaciones. El flujo magnético asociado puede tomar tres valores y = 0, +1/3, 
en función de la disposición particular de los defectos. Debido a esta dependencia con 
la disposición concreta de los defectos, es posible encontrar situaciones en las cuales 
la contribución total sea nula, no apareciendo modificación alguna en las funciones 
de onda, lo cual resulta estar en contra con simulaciones realizadas al margen de este 
modelo. 

En el capítulo cuarto presentaremos un modelo que generaliza y modifica con- 
venientemente el modelo de campos gauge analizado en el capítulo anterior. Este 
modelo se basa en la apreciación de que realmente, los electrones se hayan en grafeno 
curvado. El modelo de espinores de Dirac en un espacio curvo fue aplicado con éxito 
en el estudio de la estructura de niveles electrónicos en la molécula de Cgy y sus deri- 
vados. El fullereno, o molécula Cépy, puede considerarse como una lámina de grafeno 
en el que se han introducido cinco defectos pentagonales, formándose una estructura 
esférica. La conjunción del modelo de fermiones en una esfera con la presencia del 
campo gauge que mezcla los puntos de Dirac da cuenta de forma muy aproximada de 
esta estructura de niveles. Estando la estructura de campos gauge bien entendida, 
en esta tesis analizamos de forma profunda el efecto que tiene la curvatura sobre la 
estructura electrónica. Para ello, resumimos la teoría de espinores en espacios cur- 
vos, introduciendo como herramientas las díadas y la conexión de espín. Las díadas 
constituyen un conjunto de vectores que forman un sistema de referencia localmente 
plano en cada punto de la superficie curvada. En tal sistema de referencia es posible 
construir representaciones espinoriales de las transformaciones de simetría tales como 
rotaciones o traslaciones. Debido a que este sistema de referencia es distinto en cada 
punto, es necesario modificar el concepto de derivada para poder definir correctamente 
una acción global. 





El trabajo consiste en encontrar una métrica que describa la superficie curva, y 
tal métrica resulta ser bien conocida en cosmología, ya que los defectos topológicos 
son el análogo bidimensional a las cuerdas cósmicas en cuatro dimensiones. Estudia- 
remos configuraciones de defectos en las que de forma efectiva los campos gauge no 
geométricos (el primero y tercero de ellos explicados en el capítulo anterior). Veremos 
como podemos transformar el hamiltoniano de Dirac en espacio curvo en un hamil- 
toniano en el espacio plano en presencia de dos tipos de desorden: un campo gauge, 
que generaliza el campo gauge asociado a la matriz 2, y, lo que es más interesante, 
una velocidad de Fermi dependiente de la posición generada por el desorden, término 
que no tiene análogo en el modelo de campos gauge y que, en cierto sentido, cura las 
limitaciones de este último modelo. 

Hemos calculado perturbativamente la densidad de estados para varias disposicio- 
nes de defectos, encontrando modificaciones apreciables que pueden ser observadas 
experimentalmente, En concreto, en el caso de un defecto de tipo Stone-Wales, en 
el que el modelo de campos gauge no da alteración alguna a la densidad de esta- 
dos, se aprecia una densidad de estados que puede ser comparada con simulaciones 
publicadas en trabajos anteriores. 

Hemos mencionado que las muestras de grafeno curvadas, pueden estarlo tanto por 
la presencia de defectos topológicos o por las fluctuaciones térmicas inherentes a la fase 
plana de la lámina. Un tercer mecanismo que ocurre en las muestras depositadas sobre 
un sustrato. En tales situaciones, corrugaciones que no se encuentran suprimidas y 
las propias imperfecciones del sustrato pueden hacer que la muestra de grafeno se 
encuentre curvada de manera suave. Hemos aplicado el formalismo de espinores en 
espacios curvos a esta situación de curvatura suave, modelizada por una superficie con 
forma de gaussiana. Como en el caso anterior, encontramos modificaciones locales 
en la densidad de estados, con la salvedad de que la validez de tales modificaciones 
no se limitan a regiones concretas, al contrario que con los defectos, en donde las 
modificaciones perturbativas densidad de estados dejan de tener sentido. 

La pregunta que motiva el capítulo quinto surge de forma natural. ¿Cuál es el 
efecto, no de un conjunto pequeño de defectos, sino de un conjunto muy grande de 
ellos? Este efecto se debe apreciar, no sobre magnitudes locales, como la densidad de 
estados local, sino sobre propiedades medidas a escalas macroscópicas, tales como la 
conductividad eléctrica. Antes de calcular esta magnitud, definiremos con precisión 
el modelo y la manera de calcular observables promediados sobre todas las posiciones 
de desorden. 

En el capítulo anterior, vimos que el desorden topológico es un desorden de largo 
alcance. Esta cualidad se va a transmitir en el modelo de desorden, apareciendo en 
forma de una interacción de cuatro fermiones de tipo corriente-corriente. el parámetro 
de acoplo adimensional de dicha interacción g(k, k”) resulta ser divergente en el límite 
k' > k. Otra característica de este modelo, que resultará ser clave, es la aniso- 
tropía de la función g(k,k'). Descompondremos la interacción de cuatro fermiones 
en un conjunto infinito de canales angulares, de tal manera pasaremos de tener una 
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función g dependiente de los ángulos de dispersión de entrada y de salida a un con- 
junto de interacciones etiquetadas por un momento angular. La divergencia en el 
parámetro de acoplo hará que las propiedades a una partícula, tales como la vida 
media 1/27 o la densidad de estados al nivel de Fermi p(w = 0) sean magnitudes 
divergentes. Daremos cuenta de tal divergencia introduciendo un parámetro adimen- 
sional 0, cuyo significado físico estaría relacionado con una distancia finita sobre la 
cual, los electrones serían insensibles a los defectos. De esta manera, las propiedades 
citadas anteriormente pasan a ser finitas pero dependientes de este parámetro, que 
es esencialmente incontrolable. Emplearemos como herramienta de trabajo la teoría 
de los modelos sigma no lineales. En esta formulación, se transforma el modelo de 
fermiones de Dirac en interacción a un modelo no lineal de matrices, cuyo significado 
físico es el de der los grados de libertad de difusión del sistema. En el caso concreto de 
una interacción de tipo corriente-corriente como la que aquí nos atañe, veremos que 
aparecen dos conjuntos de modos de difusión, que a orden arbol, no interaccionarán 
entre sí, generando dos posibles canales para la difusión. 

Sin embrago, bajo las circunstancias de una densidad de estados y una vida 
media finitas, calcularemos la conductividad DC,opc, a nivel semiclásico (Aproxi- 
mación de Born Autoconsistente) encontrando que dicha magnitud es independiente 
del parámetro 0, estando pues bien definida. El motivo principal por el cual obtene- 
mos una conductividad bien definida es la existencia de múltiples modos angulares 
para la interacción, que generarán acoplos extra en las acciones para los modos de di- 
fusión. Integrando los modos de difusión masivos, encontraremos que la constante de 
difusión se ve renormalizada de tal manera que podemos obtener una conductividad 
independiente del parámetro introducido para cancelar la divergencia. 

Junto con la independencia con el parámetro 0, y lo que es más interesante, en- 
contramos que la conductividad DC es dependiente de la densidad de desorden. Esta 
dependencia, que es propia de los sistemas difusivos, choca con gran parte de los 
resultados teóricos y experimentales encontrados en la literatura hasta ahora, los 
cuales evidencian un valor para la conductividad DC universal, similar al valor que 
presentaría el grafeno en ausencia total de desorden. 

Al final del capítulo comentaremos los rangos de validez de la teoría semiclásica, 
rangos que se encuentran definidos en términos del parámetro 0. También comenta- 
remos las posibles direcciones futuras de este trabajo. 

El capítulo sexto no guarda relación directa con el trabajo anterior, versa sobre la 
posibilidad de obtener haces polarizados en el grado de libertad de valle empleando 
la distorsión trigonal que aparece a energías moderadas en las bandas T del grafeno. 

La “valletrónica” es un término que acuñaron los autores de la referencia [1] al 
proponer que el grado de libertad de valle en grafeno (es decir, la estructura de 
doble cono en la relación de dispersión a baja energía) podía ser utilizada de la 
misma manera que el espín real para transportar información, tal y como ocurre 
en espintrónica. En general, la función de onda a baja energía para un electrón 
Tr en grafeno es una combinación lineal de las funciones de ondas en ambos conos, 





apareciendo esta degeneración como un factor multiplicativo en la mayoría de las 
circunstancias. Sin embargo, es posible separar electrones procedentes de distintos 
conos o valles cuando se consideran geometrías especiales, tales como un contacto 
cuántico inmerso en una banda de espesor finito de grafeno. Esta geometría tiene el 
principal inconveniente de ser experimentalmente irrealizable a día de hoy, ya que este 
dispositivo necesita un tipo concreto de bordes (bordes zig-zag) que puedan contener 
modos cero, modos de energía cero responsables de generar un desequilibrio en la 
población de ambos conos y por tanto, de una corriente polarizada en el valle. 

Nuestra propuesta es utilizar como dispositivo una unión n —p=—n”. Esta unión, 
fácilmente obtenible experimentalmente, consiste simplemente en una barrera de po- 
tencial. La novedad consiste en considerar en uno de los lados de la barrera un 
potencial tal que la distorsión trigonal de la banda sea apreciable. Esta distorsión, 
correspondiente a considerar el siguiente término en el desarrollo en serie de momen- 
tos de las bandas, hace que las superficies de energía constante, que en el caso de 
Dirac eran circunferencias perfectas, ahora posean forma triangular, rompiéndose la 
isotropía del medio, y lo que es más interesante, la anisotropía en cada cono es di- 
ferente. La transmisión de electrones a través del grafeno en estas condiciones de 
anisotropía selectiva entre valles hace que la corriente transmitida tome direcciones 
diferentes para cada valle, produciéndose un fenómeno de separación de haces pola- 
rizados en el valle, siendo la separación angular de ambos haces, en la circunstancia 
más favorable, de aproximadamente 20 grados. 

El efecto de separación de haces no es el único que encontramos, ya que depende de 
la orientación relativa entre la barrera de potencial y la red. El otro efecto extremo, 
también consecuencia de la distinta forma de las superficies de energía constante, 
consiste en una colimación hacia delante de uno de los dos conos, y una dispersión 
angular del haz correspondiente al otro cono. Aunque cualitativamente la respuesta 
de cada cono en este caso sea muy diferente, al estar ligada a la curvatura de la 
superficie de Fermi, la diferencia de energías necesaria entre una parte de la barrera 
y la otra hace que la curvatura de la superficie de Fermi en la zona de distorsión 
trigonal sea lo suficientemente pequeña como para que esta diferencia de efectos no 
sea muy pronunciada. 

En el capítulo primero mencionamos que la aproximación de Dirac empezaba a 
fallar y el efecto de distorsión trigonal empezaba a ser apreciable en torno a energías 
de 0.6eV. Los dispositivos experimentales disponibles hoy en día no son capaces de ir 
mucho más allá de esas energías. El principal motivo por el cual no se pueden obtener 
dopados a mayor energía es debido al espesor concreto del sustrato de 510, ya que su 
valor viene marcado por su contraste óptico con la lámina de grafeno, lo que hace que 
sea observable a través de un microscopio óptico. Existe, sin embrago, otro máximo en 
el contraste óptico correspondiente a un valor mucho menor del espesor del sustrato, 
lo cual permite alcanzar valores iguales y superiores al valor citado anteriormente 
para la energía, siendo posible entonces observar experimentalmente este fenómeno 
de polarización de valle. 
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Este fenómeno de polarización de valle, aunque espectacular, no es exclusivo del 
erafeno. Ciertos cristales fotónicos (estructuras periódicas de materiales con distinto 
índice de refracción) tienen en su estructura de bandas conos similares a los conos 
de Dirac, e igualmente, zonas de energía mayor análogas a las zonas de distorsión 
trigonal en las bandas del grafeno. Al final del capítulo presentaremos simulaciones 
numéricas en cristales fotónicos de los dos efectos descritos anteriormente, simulando 
la unión p — n7” por medio de la unión de dos materiales con índices dde refracción 
diferentes. La ventaja actual en los cristales fotónicos es que este tipo de dispositivos 
experimentales ya son de uso común en los laboratorios de óptica. 


Capítulo 2 


Grafeno: un material 
verdaderamente bidimensional. 





Figura 2.1: Imagen de una muestra de grafeno tomada mediante la técnica de mi- 
croscopía de fuerzas atómicas. La lámina de grafeno es de un átomo de espesor y de 
aproximadamente 10 micras de longitud en ambas direcciones, conteniendo la muestra 
del orden de billones de átomos de carbono. 


2.1 Las formas alótropas del carbono. 


El átomo de carbono tiene una estructura electrónica 15?25?2p?. Esta estructura 
electrónica confiere al carbono la capacidad de conectarse entre sí de múltiples mane- 
ras y formar una variedad considerable de sólidos, tanto moleculares como cristalinos, 
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por medio del proceso de hibridación.Los tres tipos principales de hibridación del car- 
bono dan lugar a materiales diferentes: si en la hibridación participa un electrón en 
el orbital 25 y un electrón en el orbital 2p, la hibridación es de tipo sp!, y se ge- 
neran estructuras moleculares cuyo miembro más simple es el acetileno (C2H3). El 
papel que desempeñan cada uno de los orbitales en esta hibridación es diferente: el 
orbital hibridado cd es el que confiere estabilidad estructural a la molécula uniendo 
los átomos de carbono, mientras que el resto de orbitales p, son perpendiculares al 
orbital d se encuentran débilmente ligados. Las cadenas generadas por los átomos 
de carbono unidos por orbitales d pueden ser tan largas como se desee, haciendo 
que el número de moléculas unidimensionales sea virtualmente infinito. La siguiente 
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Figura 2.2: Los diferentes tipos de hibridación en el átomo de carbono: (a) Hibri- 
dación sp?. (b) y (c):Hibridaciones sp? 

y sp”. 
estructura posible para el carbono procede del tipo de hibridación sp?. En este tipo 


de hibridación, el orbital atómico 2s se mezcla con dos orbitales 2p para generar tres 
orbitales moleculares d, dispuestos en forma trigonal. La estructura generada es una 
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red de tipo hexagonal, conocida como grafeno. El orbital 7 restante es perpendicular 
a los orbitales d, y, al igual que en el caso del acetileno, el electrón que se encuentra 
en dicho orbital está débilmente ligado y puede saltar de orbital a orbital, y será el 
responsable de las propiedades electrónicas de baja energía de este sistema. 

Por último, la estructura restante procede de la hibridación del orbital 2s con los 
tres orbitales 2p ( hibridación de tipo sp*). Se generan cuatro orbitales d- dispuestos 
en estructura con forma de tetraedro. El sólido basado en el carbono con esta estruc- 
tura es el diamante. Al contrario que las cadenas unidimensionales de tipo acetileno 
y la estructura bidimensional del grafeno, el diamante posee una estructura tridimen- 
sional, y puesto que no existen orbitales 7 deslocalizados, el diamante será (y es) un 
magnífico aislante. 

El grafeno, a su vez, es la estructura generatriz de otros compuestos. Por ejemplo, 
la inclusión de 12 pentágonos en la red de grafeno produce una molécula cerrada, el 
fullereno o Cgp, producido experimentalmente en 1985 por R. Curl, H. Kroto y R. 
Smalley. Otras estructuras derivadas del grafeno son los nanotubos, descubiertos en 
1952 por Radushkevich y Lukyanovich, y más tarde, en 1991, redescubiertos para 
la ciencia occidental por $. lijima,. Un nanotubo no es más que una lámina finita 
de grafeno enrollada sobre sí misma formando un cilindro. La libertad a la hora de 
enrollar la lámina de grafeno hace que existan toda una plétora de estas estructuras. 
Sin embargo, debido a la peculiar estructura de baja energía del grafeno, que veremos 
más adelante, y a las simetrías de la red hexagonal, los nanotubos puede pertene- 
cer a dos grupos: nanotubos metálicos y semiconductores. En teoría, los nanotubos 
metálicos pueden soportar corrientes eléctricas muy superiores a las soportadas por 
metales convencionales como el cobre o la plata, lo que los hace muy interesantes para 
la posible fabricación de dispositivos. Por último, el grafito, que es la estructura 





Figura 2.3: Imagen de un nanotubo. 
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Figura 2.4: Imagen de la estructura atómica del diamante. 


derivada del grafeno conocida desde hace más tiempo, está formado por una multitud 
de láminas de grafeno depositadas en forma de estratos unas sobre otras, y unidas en- 
tre sí por fuerzas de Van der Walls. No disponemos de espacio suficiente en esta tesis 
para enumerar y comentar todas las propiedades térmicas, mecánicas y eléctricas de 
los distintos tipos de alótropos del grafeno, la gran mayoría documentadas y enten- 
didas desde hace años, así como de todos los usos, tanto reales como potenciales que 
tuvieren. Es curioso cómo el grafeno, siendo la estructura madre de tantas formas de 
carbono, ha sido la última en sintetizarse. El motivo principal por el cual el grafeno 
tardó en obtenerse experimentalmente es que simplemente no se creía que pudiera 
existir en la realidad como una estructura bidimensional. 

El argumento fundamental que reside detrás de tal creencia se debe a L. Landau 
y R. Peierls, que demostraron que en un cristal estrictamente bidimensional, el valor 
medio de las fluctuaciones térmicas de los fonones escalan logarítmicamente con el 
tamaño de la muestra, divergiendo en el límite termodinámico, lo cual implica que 
el cristal bidimensional simplemente se desintegraría. Más tarde, Mermin, Wagner y 
Hohenberg extenderían este resultado a sistemas bidimensionales más generales esta- 
bleciendo que no existe orden de largo alcance en dos dimensiones. No obstante, desde 
hace casi veinte años se sabe teóricamente que, en realidad, los cristales bidimensio- 
nales pueden existir en la realidad. El motivo, que veremos de forma detallada en el 
capítulo siguiente, es que las oscilaciones en la tercera dimensión pueden acoplarse a 
las oscilaciones en la superficie definida por el cristal, aniquilando las fluctuaciones 
divergentes y, por tanto, estabilizando el cristal bidimensional. Pese a todo, el cristal 
no se encontrará en una situación completamente plana, apareciendo corrugaciones y 
defectos, tal y como se observa experimentalmente hoy en día. 

A título informativo, comentemos brevemente el método que se empleó originial- 
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Figura 2.5: Imagen de la estructura estratificada en el grafito. 


mente para aislar y caracterizar láminas de grafeno. El método de obtención no es 
realmente complicado, y consiste en depositar y separar un cristal de grafito sobre un 
sustrato hasta obtener capas cada vez más finas de grafito, apareciendo en algunos 
casos capas de grafito de espesor de un átomo. El principal problema consiste ahora 
en la identificación de monocapas[2] frente a otros fragmentos de grafito de mucho 
mayor espesor. Las técnicas de imagen como la microscopía de efecto túnel no dan 
información sobre el número de capas de las muestras. El elemento clave para la 
identificación de monocapas es que el grafeno de espesor monoatómico resulta ser 
visible cuando se deposita sobre un sustrato de Si02 de un espesor determinado][3]. 
Incluso, esta técnica resulta ser muy sensible, ya que una variación del 5% en el espe- 
sor de la parte de sustrato de 5103 haría que la muestra fuera invisible. En el último 
capítulo de esta tesis, encontraremos que, con objeto de poder acceder a energías 
donde el término de distorsión trigonal es relevante, es necesario variar el espesor de 
dicho sustrato de dióxido de silicio. Afortunadamente, el máximo en el contraste ne- 
cesario para detectar monocapas de grafeno no es único, pudiéndose jugar con varios 
espesores. Además, con el tiempo han aparecido otras técnicas complementarias a la 
microscopía óptica para determinar el número de capas en una muestra de grafeno, 
tales como la microscopía Raman o el efecto Hall cuántico. 


Curiosamente, el grafeno no es el único cristal bidimensional que se ha obtenido 
experimentalmente[4]. Empleando el método de exfoliación de cristales, se han obte- 
nido cristales bidimensionales de nitruro de boro (BN), seleniuro de niobio (NbSez), 
sulfuro de molibdeno (MoS3) y monocapas de superconductores de alta temperatura 
(BizSr2CaCuz0O,) abriendo la puerta al posible estudio experimental y teórico de 
sistemas bidimensionales. 
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El grafeno es un material estrella hoy en día en Materia Condensada, gracias a 
sus múltiples facetas. Por un lado, representa un material cristalino bidimensional, 
lo cual lo hace muy atractivo para la comunidad de física estadística, al compartir 
teóricamente muchas propiedades estadísticas con las membranas bidimensionales[5, 
6]. Como veremos a continuación, la estructura electrónica de baja energía del grafeno 
viene descrita por la ecuación de Dirac sin masa, lo que abre las puertas a todo un 
campo de estudio consistente en la implementación de la Electrodinámica Cuántica 
en (2+1) dimensiones en un sistema de materia condensada][7, 8, 9, 10, 11]. Fonones 
y espectroscopía Raman|12, 13] y efecto Hall cuántico anómalo en grafeno y la bicapa 
de grafeno|14, 15, 16] son otros de los efectos interesantes, distintos a los estudiados 
hasta ahora en semiconductores convencionales. Estos fenómenos físicos nuevos y 
fascinantes han atraído la atención de los investigadores del campo. Mención aparte 
merecen las propiedades anómalas de transporte que presenta este sistema, que se 
estudiarán en el capítulo quinto de esta tesis. 


2.2 Estructura de la banda 7 en el grafeno. Aproxi- 
mación de enlace apretado. 


Como hemos dicho antes, la hibridación sp? en los átomos de carbono genera una 
estructura bidimensional hexagonal, tal y como se ve en la figura (2.6). La estruc- 





Figura 2.6: Imagen de la estructura hexagonal, característica del grafeno. Los vectores 
de la red directa son a, = V3a (3, 3) y az = vV3a (3, 4). 


tura cristalina hexagonal puede entenderse como una superposición de dos subredes 
triangulares interconectadas, esto hace que la red hexagonal no sea una auténtica red 
de Bravais, y que la celda unidad en el espacio real esté constituida por dos átomos 
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pertenecientes a distintas subredes, separados una distancia a = 1.424. En la figura 
(2.6) se muestran los vectores de traslación a, y az para la red de grafeno, mientras 
que los vectores 0 2 3 representan los vectores a los primeros vecinos. Los primeros 
vecinos del átomo de la otra subred son simplemente —01 23. 

Vamos a calcular la estructura de bandas para los orbitales 7 en aproximación 
de enlace apretado[17], a primeros vecinos, para una lámina de grafeno perfecta e 
infinita. En tal aproximación, un electrón situado en una celda unidad viene descrita 
por la siguiente función de onda, en virtud del teorema de Bloch: 


AB N AB 
v(K FEE Ce in pl(e— Rs.) = Y CndBn(K,r). (2.1) 
m Ím m 


En la función de onda (2.1), el índice m = A, B representa cada una de las dos 
subredes, el índice j,, recorre todas las posiciones R;,, de los átomos de cada subred 
m en todo el cristal y H1p(r) es la función de ondas del orbital p en cada átomo 
normalizada. El vector K se encuentra cuantizado típicamente en unidades de E 
Conocido el ordenamiento de los átomos de carbono en la red hexagonal y las itdones 
de onda atómicas para los orbitales p, el problema consiste en calcular los coeficientes 


Cin de tal forma que la función de onda (2.1) sea un autoestado del hamiltoniano 


H(r) = Y + V(r), (2.2) 





con la propiedad de periodicidad de la red del potencial V, V(r) =V(r —R)): 
H(DV(K, r) = E(k)v(K, r). (2.3) 


Multiplicando la ecuación (2.3) por W*(K,r) e integrando en r, tenemos que 


N AB dE 

yy L 00, ¿KR Rin 2 ro —Ri,)H (Or Ri) = 

o a AB 

A OO ein) | rg (Rió Rin). (24) 
jmyin N ¿»Mm 


En este punto, consideraremos la aproximación basada en suponer que la integral en 
el lado derecho de (2.4) es no nula únicamente cuando ¿y = jm (solapamiento débil 
entre estados en distintas subredes), de tal forma que podamos escribir 

ENAB A,B 

a 7 Kimi CIC, ¿nm = E Y C%, Ci. (2.5) 


a n,m 


Más adelante comentaremos las implicaciones que conlleva esta simplicación. 
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Centrémonos ahora en el lado derecho de la ecuación (2.4). Debemos notar que 
los vectores de traslación de ambas subredes se relacionan entre sí, para 1 = j, por 


R; — R; = —(092. (2.6) 


Esta relación resultará básica en lo que sigue. 
La aproximación de enlace apretado a primeros vecinos implica que la integral 


dE Bro*(r —R,)Hó(r—R,,) = ] Bro*()Hóo(r— (Rs, —Ri,)) (2.7) 


será distinta de cero en dos posibles casos, habiendo entonces dos contribuciones 
en el lado derecho de (2.4). El primer caso corresponde a la situación en la que 
R,,, = R;,, reduciéndose entonces la integral (2.7) a e = fd%ro*(r)Hóo(r). El 
término que contribuirá al lado izquierdo tendrá la forma € 7 CE 
A la segunda contribución, correspondiente a considerar primeros vecinos ,n 4 m, 
le corresponderá una fase que puede escribirse en términos de los vectores R;, de una 
misma base, usando la relación (2.6), como 
R,, —R 


E (2.8) 


in tm 


Ahora que tenemos ambas sumas en la misma base, y debido a la invariancia de las 
funciones de onda bajo traslaciones de cualquier vector de la red real, se tiene que 


AB 3 No 
Y Cn Y ini | Pro HA(r — (Rin — Rin — 62)) = 


nám jm im 
A,B 1 2 y 

A / Pro (r) H(mo(r— (a, -aj-0)).. (29) 
nám J,1 


Ahora, simplemente, es suficiente con calcular a; — a; — 9, para 1, j = 1,2, siendo el 

resultado de las posibles combinaciones para ¿ y j los tres vectores —01 23. Dado que 
la integral f d“ró*(r) H(r)p(r — 0123) no depende del vector 4 concreto, por lo que 
definiremos el parámetro de salto t como t = f d*r4*(r) H(r)jó(r — 61 23) (t - 2.8eV). 
Finalmente, el correspondiente término en el lado izquierdo de la ecuación (2.4) queda 
de la siguiente forma: 


A,B AB 
DR du CN SR Cia E 9 (2.10) 


nám j nám 


Reuniendo todas las partes de la ecuación, podemos escribir (2.4) en forma de ecuación 


matricial: 
os de ) o ) - E(k) a ) (2.11) 
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Figura 2.7: Estructura de las bandas r en el grafeno. 


Las bandas, representadas por los autoestados de la ecuación de autovalores (2.11), 
son, empleando los valores explícitos de los vectores d; 


k k 
E(k) = 0 0] 1 + 4cos (47) COS (E .) + 4 cos? (E) (2.12) 


La estructura de bandas se muestra en la figura (2.7). En la práctica, el parámetro 
€o puede eliminarse redefiniendo el origen de energías en la ecuación (2.11). La 
zona de Brillouin en el grafeno resulta tener forma hexagonal, hecho que puede an- 











ticiparse conociendo los vectores que definen la red recíproca, b¡ = == (1, v3) y 
b; = El (1, 1/3). Tal y como puede verse en la figura (2.8), las bandas de valencia 


y de conducción, representadas por las soluciones E(k) < 0 y E(k) > 0 respectiva- 
mente, se tocan en seis puntos, correspondientes a las esquinas de la zona de brillouin. 
Por periodicidad de la red recíproca, el número de puntos independientes se reduce a 
dos: K; = : (0, 15) y K_; = : (0, 25). Este par de puntos, en los cuales los cuales 
las bandas se tocan se conocen como puntos de Dirac. Puesto que cada átomo de 
carbono cede un electrón del orbital 2p a las bandas, tendremos que la banda de va- 
lencia se encuentra totalmente llena y la de conducción totalmente vacía, situándose 
el nivel de Fermi, para el material neutro, justamente en los puntos de Dirac, y las 
excitaciones electrónicas a baja energía se encontrarán en el entorno de dichos puntos. 
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Figura 2.8: Estructura de las bandas r en el grafeno. La zona de Brillouin exagonal 
se muestra superpuesta. 


2.3 Modelo continuo de baja energía. 


Estos estados de baja energía, responsables de toda la física que se describirá en 
esta tesis, son, por lo tanto, aquellos que poseen un vector de ondas muy próximo 
a alguno de estos puntos, por lo que el hamiltoniano de baja energía que describa 
dichos estados puede obtenerse desarrollando en torno a los puntos de Dirac K; _1 
en la ecuación (2.11). Puesto que tenemos dos puntos de Dirac en torno a los cuales 
podemos desarrollar en serie, el resultado va a ser que, de forma efectiva, tenemos dos 
tipos de estados de baja energía, dependiendo de entorno a cual de los dos puntos de 
Dirac nos encontremos: 


1. a de ) li did ) +HO(k), (2.13) 


en donde el parámetro s +1 nos informa del punto de Dirac en torno al cual desa- 
rrollamos, y Up = da es un parámetro con dimensiones de velocidad, que tiene un 
valor muy aproximado de vp  10%m/s, y recibe el nombre de velocidad de Fermi. 
Podemos escribir la ecuación (2.13) de forma más compacta definiendo las matrices 


cuadridimensionales “y a partir de las matrices de Pauli ur y 7: 


y =x%r080", =1807, (2.14) 
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quedando finalmente el hamiltoniano 
H =vpyk, (2.15) 


que resulta ser formalmente el hamiltoniano de Dirac en dos dimensiones espaciales. 
Llegamos pues a la siguiente conclusión: los estados electrónicos de baja energía en 
las bandas r se pueden entender como estados que satisfacen la ecuación de Dirac sin 
masa, actuando la estructura de subred como un grado de libertad de espín. En la 
literatura sobre grafeno, es muy común llamar a dicho grado de libertad pseudoespín 
o espín AB, en referencia a la estructura de subredes de la red hexagonal. 

Los coeficientes C¿K pueden escribirse en forma espinorial, que, en el caso de baja 
energía donde hemos encontrado que el hamiltoniano tiene la forma del hamiltoniano 
de Dirac, pueden escribirse como 


¿os 
CK ect 
CK Ae? 
sl ce |- dez ll (2.16) 
Ay en 
Cr —le ci? 


El parámetro A = +1 indica a qué banda nos estamos refiriendo, A= 1 para la banda 
de conducción y A= —1 para la banda de valencia, y el ángulo 0, es el ángulo del 
vector de onda k = k (cos Ox, sin 0x.). 

Pasemos ahora a comentar las aproximaciones que hemos hecho en el modelo, que 
son básicamente dos: la primera está relacionada con el número de vecinos a los que 
permitimos que los electrones “salten”. Si incluimos por ejemplo, términos de salto a 
segundos vecinos, encontraremos un término en la ecuación 2.11 que conecta átomos 
de la misma subred, por el simple hecho de que todos los segundos vecinos de un 
átomo pertenecen a la misma subred que éste, y que tendrá estructura diagonal. Por 
otro lado, puesto que cada segundo vecino posee otro situado en su posición inversa, 
el desarrollo a bajos momentos de este término contendrá un término constante, 
que romperá la simetría entre la bandas de valencia y conducción (simetría electrón 
hueco) e incluirá un término cuadrático en k, Esy = A El parámetro t' es 
la probabilidad de salto por unidad de tiempo a los segundos vecinos, y se estima 
que tiene un valor de t' => 0.1eV. La inclusión de más vecinos en la aproximación de 
enlace apretado es posible, pero carece de interés práctico, puesto que a baja energía 
no aparecerán términos distintos a los que puedan aparecer al considerar primeros 
o segundos vecinos, y los términos que aparezcan irán acompañados de coeficientes 
análogos a t o t', pero cuyo valor será totalmente despreciable frente a estos dos. 

La segunda aproximación es la aproximación de solapamiento débil hecha en la 
ecuación (2.5), donde consideramos que la integral 


Jero (r—Ri)ó(r Ry) 
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era nula salvo cuando R,, =R,,,, es decir, en el propio átomo en cuestión. Podemos 
perfectamente considerar la posibilidad en la que dicha integral sea distinta de cero 
cuando los orbitales pertenezcan a átomos adyacentes, o más allá. La aproximación 
de ser no nula a átomos adyacentes conlleva modificar el lado derecho de la ecuación 


(2.11) por el término matricial 


Elk) A ) | 


tal y como se desprende sencillamente de la ecuación (2.4). El solape entre orbitales 
de átomos vecinos es a, y se ha estimado que tiene un valor en torno a 0.13. La 
relación de dispersión resultante, 


_ cost (K) 
LEO" 


como en el caso anterior, rompe la simetría electrón-hueco. También, como en el 
caso anterior, las modificaciones que conlleva esta aproximación solamente se ven a 
segundo orden en el desarrollo en serie de f(k), siendo el término lineal en el vector 
de onda el mismo en todos los casos, y el descrito en la ecuación (2.15). 

Como vemos, la ecuación de Dirac (2.15) resulta ser una ecuación de ondas de 
baja energía muy robusta frente a las aproximaciones que se hicieron para obtenerla 
(ya que, aunque no lo hayamos comentado antes, la forma de (2.15) está protegida 
por las simetrías de la red cristalina). Es posible ver en experimentos de fotoemisión 
resuelta en ángulo[18] en qué rango de energías es válida esta aproximación de Dirac. 
A simple vista, podemos estimar que la aproximación lineal en las bandas es válida 
hasta una energía del orden de 0.6 —0.7eV, energías en las que el término cuadrático, 
también llamado de distorsión trigonal, comienza a ser relevante. A lo largo de esta 
tesis, exceptuando el último capítulo, donde estudiaremos un efecto de la distorsión 
trigonal, vamos a considerar exclusivamente la aproximación de Dirac. El valor de 
la energía nos introduce una escala de distancias mínima para la cual la ecuación 
de Dirac describe la física de forma precisa. Escribiendo la velocidad de Fermi en 
unidades de Ármstrong y electron voltios ( up — 6.5AeV), dicho valor mínimo es del 
orden de 114, es decir, la ecuación de Dirac va a describir perfectamente los estados 
electrónicos en escalas que puedan ir desde miles de Ármstrongs a 10 = 11 Á, siendo 
entonces una magnífica aproximación para describir todos los efectos macroscópicos 
que deseemos estudiar. 

En los capítulos siguientes vamos a estudiar la densidad de estados del grafeno en 
presencia de desorden topológico, por lo que es conveniente conocer la densidad de 
estados en el grafeno perfecto. 

Se define la densidad de estados como el número de estados electrónicos por unidad 
de energía y unidad de volumen|19): 


pl / Pkó(w — E(k)). (2.17) 


E(k) 
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Figura 2.9: Densidad de estados en las bandas 7 de conducción y de valencia en 
grafeno. La densidad de estados en el nivel de Fermi w = 0 es nula. 


En la figura (2.9) mostramos la densidad de estados para la relación de dispersión 
(2.12). El elemento más relevante a destacar de la densidad de estados en las bandas 
TT es que es nula en el nivel de Fermi. Esta propiedad va a ser responsable o al 
menos estará involucrada en todas las propiedades anómalas que exhibe el grafeno, 
desde el punto de vista de la física de semiconductores. La aproximación de Dirac 
nos va a capturar correctamente la densidad de estados en torno al nivel de Fermi. 
Efectivamente, utilizando la aproximación E(k) = +up|k| en (2.17), se encuentra 
fácilmente que la densidad de estados en la aproximación de Dirac es lineal con la 
energía: ; 

pole) = la (2.18) 

En el lenguaje de segunda cuantización(20], la densidad de estados local es la parte 

imaginaria de la función de Green a dos puntos 


plw,r) = lim ImTrG(w, r,r”). (2.19) 


Con la definición de la densidad de estados, finalizamos la descripción breve del 
grafeno y de su estructura electrónica a baja energía. En los siguientes capítulos, nos 
centraremos en el estudio de estas propiedades en presencia de desorden topológico. 


Capítulo 3 


Fermiones de Dirac en grafeno con 
defectos topológicos. 


“Like human defects, those of crystals come in a seemingly endless variety, 
many dreary and depressing, and a few fascinating. [19/” 


3.1 Introducción. Defectos en membranas bidi- 
mensionales. 


En la vida real, los cristales tienen defectos que están presentes incluso en equilibrio 
térmico. En muestras reales de grafeno pueden encontrarse todo tipo de defectos: 
vacantes, defectos intersticiales, microroturas, átomos de otras especies, impurezas 
magnéticas y dislocaciones. Es bien conocida la importancia que tienen los defectos 
a la hora de entender las propiedades de transporte electrónico y las propiedades 
mecánicas de un material. En el caso del grafeno, el papel que puede jugar los 
defectos topológicos es incluso más interesante, al ser un material auténticamente 
bidimensional. 

A continuación, revisaremos de forma breve el problema de una membrana cris- 
talina bidimensional. Veremos que si se permite a esta membrana tanto deformarse 
como combarse, a temperaturas suficientemente bajas la membrana es capaz de esta- 
bilizarse en su fase plana[21, 22, 23]. 

Supongamos que la posición sobre la superficie de la membrana viene descrita por 
el vector tridimensional u(x,y) = (1,y,0) + (u(x, y), uz(x, y), f(x, y)) como función 
de las coordenadas del plano (x, y), por lo que las distancias sobre la superficie dis- 
torsionada verifican 

lr == Bro + ZUALALE (3.1) 


El tensor u,¿ está definido, en orden más bajo en gradientes de u; y f, como 
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Nótese el término anarmónico en (3.2), que acopla los modos 0;u, con 0;f. 
En la teoría de la elasticidad, la energía necesaria para deformar la membrana es 


Elu, f] = > Jal 2quuz, + Mia) + 25 fal (vay. (3.3) 


Las constantes y Á son las constantes de Lamé de la membrana, y K es la rigidez 
de doblado del material. El parámetro de orden que vamos a emplear para definir el 
orden cristalino es el vector normal a la membrana, n(r), relacionado con la función 


f a través de 
1 


ETA 


El valor medio del vector n puede calcularse a partir de (3.3) empleando toda la 
maquinaria de física estadística: 


n(r) = -0f,-0,f, 1) (3.4) 


E 


DuDfn(r)e FT. (3.5) 


La ausencia de orden cristalino se da, como sabemos de física estadística, cuando las 
fluctuaciones en torno al valor medio (3.5) son lo suficientemente grandes, haciendo 
que la fase plana de la membrana sea inestable. Las fluctuaciones en torno a (3.5) 
pueden calcularse por medio del ángulo O que forma n con el vector normal al plano 
2124): (8%) = (IV). 
Puesto que n(r) es un funcional únicamente de la función f, podemos integrar 
los modos u;¿ en la función de partición definida a partir de (3.3), obteniendo una 
energía libre efectiva 


00; 
Ej=58 [avr + K fa (3 (005 FDO n). (3.6) 
siendo K el módulo de Young bidimensional. El término no lineal en (3.6) hace que 
la rigidez de doblado deje de ser constante, pudiendo ser calculada como una serie 
de potencias del módulo de Young K. Pasando al espacio de momentos, es posible 
ecribir una ecuación de autocosistencia para la rigidez de doblado renormalizada 





dk qi(0;; — e Ja; 
47? kp(k)|q ma k]+" 


cuya solución es, en el límite de pequeños valores de q, 


«r(q) = ykgTKq 7. (3.8) 


Las fluctuaciones (0?) tienen entonces la forma 


kgT y d%q1 
2 PO pe q 
a = tur [E 412 kg(q Sí > N K ) 4rq oe 





ar(q) =«+kgTK ] (3.7) 
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La ecuación (3.9) es divergente únicamente en el ultravioleta, siendo finita en el in- 
frarrojo. Puesto que estamos trabajando con membranas cristalinas, la divergencia 
ultravioleta se suprime imponiendo un valor máximo para el momento q del orden 
del inverso del espaciado de la red, A = z Que las fluctuaciones (3.9) no sean di- 
vergentes infrarrojas significa que no son capaces de destruir la fase plana, al menos 
a temperaturas bajas (a temperaturas suficientemente altas, la entropía de la mem- 
brana cristalina favorece que ésta se arrugue, lo que sugiere una transición de fase 
desde la fase plana a la fase arrugada). El resumen de la discusión de los párrafos 
anteriores es que la combinación de los modos de deformación en la superficie,u;;, 
con los modos fuera del plano o flexurales, f, hacen que la fase plana se estabilice a 
bajas temperaturas, permitiendo pequeñas fluctuaciones de la membrana, pero man- 
teniéndose el orden cristalino. Por lo tanto, haciendo extensiva la discusión anterior 
al grafeno, podemos pensar que una lámina de grafeno suspendida tiene como fase 
estable, a bajas temperaturas, la fase plana. 


LSD» : 
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Figura 3.1: Defectos pentagonal y heptagonal en una membrana triangular. Se ha 
permitido que la membrana se curve en torno a los defectos en la tercera dimensión. 


3.2 Defectos topológicos. 


El análisis anterior nos sugiere que una muestra de grafeno libre puede mantenerse 
en su fase plana, permitiendo fluctuaciones sobre dicho estado que harían que la 
membrana se curvara localmente en regiones locales. Otra posible fuente de curvatura 
en la membrana de grafeno es la presencia de defectos topológicos: disclinaciones y 
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dislocaciones. Este tipo de defectos se al sustituir un anillo hexagonal por un polígono 
diferente. 

En las figuras (3.1) y (3.2) se muestran imágenes de discliniaciones positiva y ne- 
gativa y una dislocación, respectivamente. Una forma de construir una disclinación es 
por medio de la construcción de Volterra[25]: Supongamos una superficie cristalina 
bidimensional sobre la que retiramos o añadimos una sección angular de material. 
Uniendo los bordes obtenemos, o bien una superficie cónica, si hemos retirado mate- 





Figura 3.2: Par heptágono-pentágono curvado en la tercera dimensión. 


rial, o bien una superficie hiperbólica, en el caso de añadir material a la superficie. 
Una dislocación puede entenderse como un defecto compuesto de dos disclinaciones 
cercanas. Las disclinaciones se caracterizan por el ángulo de material retirado o 
añadido 0(r), mientras que las dislocaciones vienen caracterizadas por su vector de 
Burgers b(r), que se define como el vector que une los núcleos de las disclinaciones 
que forman la dislocación. La característica fundamental que van a cumplir las dis- 
clinaciones está relacionada con el exceso o defecto angular antes citado: cuando se 
considera un circuito cerrado C' que rodea la disclinación, el recorrido angular va a 
ser diferente al correspondiente para la lámina plana: 


pB 27 
$ 00) =1- =9 E. (3.10) 

Estos defectos juegan un papel fundamental en las transiciones de fase en sistemas 
bidimensionales[26, 27]. 

Al contrario que las ondas elásticas, que son excitaciones elementales de los cam- 
pos u¡¡ y f, y que pueden excitarse con energías arbitrariamente pequeñas, los de- 
fectos topológicos se caracterizan por tener energías de formación fijas. A partir de 
las energías de deformación (3.3) o (3.6) podemos obtener las ecuaciones de movi- 
miento para el campo f. Las ecuaciones extraídas de (3.3) son las ecuaciones de 


24 Capítulo 3: Fermiones de Dirac en grafeno con defectos topológicos. 





von Karman[28), un conjunto de ecuaciones diferenciales en derivadas parciales no 
lineales. Estas ecuaciones pueden reducirse escribiéndolas en función de la función de 
tensión de Airy x(r) y del campo f(r), teniendo que cumplir las soluciones corres- 
pondientes a los defectos topológicos la condición (3.10). Una vez encontradas dichas 
soluciones, las energías elásticas para las disclinaciones pentagonal y heptagonal son 
Es = ¿kIn(R/a) y E, = rrIn(R/a). En estas dos expresiones vemos que las energías 
divergen con el tamaño de la muestra R, haciendo que la formación de disclinaciones 
aisladas en membranas de tamaño macroscópico sea muy costosa energéticamente. 
Sin embargo, éste no va a ser el caso de las dislocaciones. 


En la figura (1.3) se muestra el valor de la energía asociada a un par heptágono- 
pentágono en unidades del módulo de Young Ky como función de la coordenada 
radial R para varios valores del cociente entre Ky y la rigidez de doblado k. Salvo 
en el límite Ky/£ —= 0, correspondiente a la situación de lámina cristalina, en la 
que no son posibles las deformaciones, y que como sabemos, hace que F" = log R, la 
energía de formación tiende a un valor finito cuando R — oo. Este hecho genera un 
cambio profundo en la estadística de formación de defectos. En el umbral en el que 
las deformaciones locales de la red en torno a los defectos son posibles, la contribución 
entrópica domina a cualquier temperatura distinta de cero. La densidad de defectos 
es ahora finita: 


np = a e FMesT) (3.11) 
De la ecuación (3.11) podemos extraer la longitud de correlación translacional £p: 


Ep np e aer /erT), (3.12) 


Esta relación nos dice que el orden traslacional estará roto hasta distancias iguales 
o mayores que £p. La ecuación (3.11) también nos dice que la densidad de defectos 
compuestos por un par heptágono-pentágono será finita ( al serlo F') a temperaturas 
finitas, en contraposición con lo que ocurre en el límite de lámina cristalina o en la 
situación de defectos aislados, en las que dicha densidad es nula en el límite R = 0. 
En el caso del grafeno, el módulo de Young toma valores típicos del orden de 1000 
GPa mientras que la rigidez de doblado k es del orden de 1 eV en nanotubos|29]. Esto 
hace que, en unidades adecuadas, da < 1, alcanzándose el valor asintótico constante 
rápidamente. 


En base a toda la discusión anterior, es muy razonable pensar que los defectos to- 
pológicos son una fuente natural de desorden en el grafeno. A lo largo de lo que queda 
de capítulo y en los siguientes, veremos la manera profunda en la que la presencia de 
este desorden modifica la estructura electrónica de baja energía y las propiedades de 
transporte electrónico de una lámina de grafeno. 
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Figura 3.3: Energía de formación de un par heptágono pentágono para diferentes 
valores de las constantes de rigidez en el caso de una membrana de estructura cristalina 
triangular. En todos los casos la energía de formación tiende a un valor constante en 
el límite infinito de radio de la muestra R salvo en el límite inextensional. El par se 
encuentra separado a una distancia de uno, en unidades de espaciado de la red. 


3.3 Efecto de una disclinación aislada. 


Podemos construir un primer modelo para describir el efecto de una disclinación 
aislada en una lámina de grafeno, basándonos en la transformación de las condiciones 
de contorno que impone el cortado y pegado de la red tras la inclusión de un defecto 
en campos gauge. Esta transformación puede considerarse, de cierta manera, como el 
procedimiento inverso al efecto Aharonov-Bohm, en el que el efecto de un campo gauge 
se transformaba en una condición sobre las funciones de onda. Este modelo se propuso 
inicialmente para explicar el espectro discreto de baja energía de los fullerenos[9], y 
posteriormente para entender el efecto de dichos defectos en nanoconos de grafito[30, 
3: 

Cuando consideramos únicamente una disclinación aislada, el efecto de cortado 
y pegado de la red hexagonal va a generar dos condiciones de contorno no triviales 
sobre las funciones de onda. Como se muestra en la figura (3.3), la inclusión de un 
defecto puede entenderse fácilmente por medio de cortado y pegado. En el caso de un 
defecto pentagonal, se ha retirado un sector correspondiente a un ángulo de P = 7/3, 
y posteriormente se han unido los átomos de los bordes. El primer efecto que genera 
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Figura 3.4: Izquierda: Efecto de un defecto pentagonal en la red de grafeno. Derecha: 
Procedimiento de corte y pegado de la red plana de grafeno para formar dicho defecto 
pentagonal. 


una condición de contorno es que ahora el recorrido angular que, desde un punto 
de vista clásico, debería hacer un electrón al rodear el defecto, es menor que en la 
situación plana, exactamente el valor angular del sector retirado. Mientras que en la 
situación plana, la función de onda satisface la siguiente condición de contorno 


V(r,0) = —V(r,0 + 27) = exp(21:2)V (r, 0), (3.13) 


siendo Y = 373 8 03 = y el generador de rotaciones en torno al eje perpendicular a 
la lámina de grafeno, en el caso que nos ocupa, la anterior condición de contorno se 
modifica de forma trivial considerando como la nueva variable angular 0” = 0(1 — 6) 


W(r, 0) = exp (2. ( — z) 2) v(r,0'). (3.14) 


Es decir, aparece una condición de contorno extra con respecto a la situación plana 
en forma de Ry = exp(-¡£Y). 

La segunda condición de contorno es ligeramente más sutil. Cuando se incluye 
un defecto cuyo número de lados es impar, como es el caso de un pentágono o un 
heptágono, los átomos en los bordes y que ahora se encuentran enlazados pertenecen 
a la misma subred, mientras que, como ya sabemos, en la situación plana cada átomo 
se encuentra unido a átomos de la subred contraria. 

Consideremos, para entender de manera más sencilla cómo describir este efecto 
como condición de contorno, la inclusión de un triángulo, tal y como se expone en 
la referencia [32], en ambas situaciones, como hemos dicho, la condición de contorno 
será la misma. En la figura (3.5) se muestra la situación y la equivalencia entre los 
átomos de ambos bordes. También se muestra como tras realizar dicha equivalencia, 
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un desplazamiento a lo largo de la costura aparece como una inversión, que genera 
el intercambio de las subredes. Esta inversión de la red tiene como consecuencias 
asociadas, en el espacio de momentos, el cambio del vector k por —k, y por tanto 
el intercambio de los puntos de Fermi, que satisfacen esa misma relación: K' = 
—K. Por tanto, si consideramos que en O = 0 el espinor tiene estructura Y = 
(Pax, Ver, Varo, Vero), al realizar una rotación de ángulo 7 y cruzar la costura, 
el espinor tendrá la forma W* = (Wgkx», VAk”, UBk, bar”. La forma matricial de la 
condición de contorno que cumple RpV = ¡VW* es 


Rp =i il y (3.15) 


El factor ¿ se incluye en la condición de contorno para que se cumpla R4 = 1,manteniendo 
la condición de contorno espinorial (3.13). 





Figura 3.5: Defecto pentagonal de apertura 6 = Tr. En la figura de la izquierda se 
muestra una trayectoria rectilínea al cruzar la costura en términos de celda unidad 
de la red hexagonal. En la figura de la derecha se muestra la superfice cónica y la 
misma trayectoria. Se observa claramente la inversión efectiva de la red con respecto 
a la trayectoria al cruzar la costura. Esta inversión de la red es la responsable del 
intercambio efectivo de los puntos de Fermi, tal y como se explica en el texto. 


Es fundamental notar que estas condiciones de contorno hacen que la función de 
onda electrónica sea discontinua en Q = 0. Puesto que en la realidad, el proceso de 
cortado y pegado que genera las condiciones de contorno no es más que un artefacto 
ideado para analizar el efecto de estos defectos sobre funciones de onda, éstas no 
tienen por qué ser discontinuas, lo que motiva la conversión de dichas condiciones de 
contorno en campos gauge, los cuales sí permiten una transformación continua de la 
función de onda a lo largo de todo el espacio. Para tal fin, el primer paso es escribir 
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las condiciones de contorno en forma de exponenciales. En el caso de Ry el trabajo 
ya está hecho en la ecuación (3.14). Para Rp tenemos|[9, 30, 33] 


o: 27 0 (05) 
Rp =1exp (+ e )) l (3.16) 


Ahora podemos entender las condiciones de contorno de forma general como 
v(2r) = Ev (9) = PAR EGO), (3.17) 


donde queda definido un campo gauge no abeliano local Ap, (1) = A(r)R; = US (r)VU;(r) 
asociado a cada condición de contorno. Por lo tanto, podemos escribir la función de 
onda general como 

V(r) = e **Dy(r), (3.18) 


con O(r) = Zp, $, Ar, (1)dr y W'(r) una nueva función de onda que satisface la 
condición de contorno periódica para el grafeno plano (3.13). 

Podemos aplicar la propia transformación gauge al hamiltoniano de Dirac para 
obtener el hamiltoniano para W'(r): 


H' =P ge? — iypy (a —i ( E z) YAy(r) — ¡MrA ml) , (3.19) 


donde, recordemos, = = 373 e 1 y Rp = —712 8 07. Podemos completar la trans- 

formación de las condiciones de contorno a los campos gauge definiendo los flujos 
.o. . e B a 1 pol 1 

ficticios asociados a cada campo gauge como Py = ¿¿ = ¿ y Pp, = 3, en caso de 

considerar un defecto pentagonal. 


3.4 Múltiples defectos. Limitaciones del modelo 
de campos gauge. 


Nuestro objetivo final es poder describir la situación en la que la muestra de 
grafeno contiene un numero arbitrario de defectos topológicos, repartidos de manera 
arbitraria por toda la muestra. La forma de extender el formalismo anterior al caso 
de múltiples defectos situados en posiciones arbitrarias es directa. Consideremos un 
punto al azar de la muestra como origen ( al final veremos que el resultado es inde- 
pendiente de esta elección). Respecto a dicho centro, las posiciones de N defectos 
serán [r,,r2,...,ry). Para ver cual es la condición de contorno no trivial inducida 
por los N defectos y, por tanto, el campo gauge ficticio asociado, consideraremos un 
circuito cerrado centrado en r = 0 que englobe a todos los defectos y a lo largo del 
cual integraremos el campo gauge ficticio total que dé como resultado la condición 
de contorno deseada. La condición de contorno no debe depender de la forma con- 
creta de dicho contorno, por lo que podemos deformarlo de tal manera que podamos 
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Figura 3.6: Deformación del contorno C.. 


aprovechar la discusión de los párrafos anteriores y entender la fase en función de las 
fases aisladas generadas por cada defecto. 

Consideremos como ejemplo la deformación del circuito C' propuesta en la figura 
(3.3). El contorno deformado consiste en traslaciones desde el punto tomado como 
origen hasta la posición de cada defecto, rotaciones en torno a cada defecto y de nuevo 
una traslación al punto de origen. El siguiente paso es ver cómo se transforman las 
funciones de onda por este camino deformado. En los párrafos anteriores hemos visto 
el caso de los desplazamientos a lo largo de caminos cerrados rodeando los defectos. 
Veamos ahora las traslaciones. El operador de traslación para la función de onda en 
aproximación de enlace fuerte es[31] 


T(n, m) = exp (Ein — m)r3 O 1) : (3.20) 


La forma de (3.20) merece un comentario. Sabemos que la función de onda total para 
un electrón en grafeno, en aproximación de enlace apretado, es 
Votar (1) == y y Tor A), (3.21) 
j=A,B i=K,K' 
en donde Vj¿(r) x eX"t, son las funciones envolventes que satisfacen la ecuación de 
Dirac, y F; (1) = Y e Y ¿(1 — d;) son las funciones de onda de Bloch. Al aplicar una 


traslación 75 de recorrido, digamos, 4 = may + naz, la función de onda Vota (1) se 
modificará con respecto a (3.21) de la siguiente manera: 


VW iotal (r + 9) = T(9)V totar (1) = eE a (1). (3:22) 
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La fase e** es global a todas las componentes del espinor, siendo por tanto, en notación 
matricial, proporcional a la matriz identidad, y, tras completar el circuito deformado 
C, su contribución será nula. Sin embargo, en el caso de la fase eXi? la situación no 
va ser tan simple, ya que, empleando los valores concretos de los puntos K y K' en la 
misma notación quiral que 4, obtenemos (3.20). Por supuesto, al final del circuito C 
la contribución de la fase como función de d debe ser nula, como en el caso anterior, 
pero vemos en (3.20) que T'(9) es proporcional a la matriz 1 € 73, y por tanto, en la 
secuencia ordenada que constituye el circuito C*: 


C=T(r,)oP,oT"*(r1)0...oT(ry)oPyoT *(ry), (3:23) 


T(r5) y Rp, no conmutarán, al depender de las matrices 73 9 1 y 07 8 72 respectiva- 
mente. El operador asociado a la curva C entonces tiene la forma general, tras operar 
adecuadamente, asociando todas las fases de la misma naturaleza: 


Po = (Rp, exp (2 Bi(N5 — N,)D) exp E le — my)73 O y] A (3.24) 
j 


expresión escrita en función del número de defectos pentagonales N5 y heptagonales 
N” presentes en la muestra, y del numero total de ellos: N = N;¿ + N7. Los dos 
primeros términos de (3.24) son básicamente las fases ya descritas anteriormente. El 
tercer término corresponde a una nueva fase, cuyo origen es, como hemos visto, la 
no conmutatividad de los operadores de traslación y de rotación en torno a cada 
defecto[31]. El flujo magnético ficticio asociado a este tercer campo puede leerse 
directamente de (3.24), dz = 30 (nm; — mj). Esta suma da lugar a tres posibles 
valores módulo 27: O y +3. 

Como discutimos en el apartado anterior, uno de los defectos más probables desde 
el punto de vista energético es el defecto de tipo Stone-Wales, consistente en dos 
parejas de heptágono-pentágono. Asimismo, en las láminas planas de grafeno se 
observan irregularidades con curvatura local que pueden deberse a la presencia de un 
número igual de defectos heptágono-pentágono. Es interesante preguntarse, llegados a 
este punto, sobe cual es el efecto de tener el mismo número de defectos pentagonales y 
heptagonales en una muestra, con el propósito de mantener la situación de planaridad 
a escala macroscópica. En la ecuación (3.24) dicha situación se alcanza haciendo de 
manera inocente N5 = N7 (y por tanto siendo N = N; + N; un número par). Esta 
simple substitución hace que los flujos asociados a los campos Ap y Ay sean nulos. Si, 
además, la configuración de los defectos es tal que el flujo Pz fuese nulo, se llega a la 
conclusión de que dicha configuración es equivalente a la situación de no tener ningún 
tipo de defectos. Esta situación es la que se daría en el caso en el que todos los defectos 
pentagonales y heptagonales estuvieran organizados en defectos Stone-Wales[34]. 

Esta situación parece estar en contradicción con estudios teóricos existentes en la 
literatura, tanto sobre láminas de grafeno|36, 37, 38] como de nanotubos[35] (figura 
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Figura 3.7: Esquema de un defecto Stone- Wales, constituido por dos pares de defectos 
pentagonales y heptagonales adyacentes, y de sus posibles orientaciones relativas a la 
red de grafeno. Estos defectos son de particular importancia en la realidad, ya poseen 
la energía de formación más baja de entre todos los posibles defectos que puedan ser 
construidos combinando defectos aislados. 


3.8), en los cuales se reporta una densidad de estados en presencia de estos defectos, 
diferente a la correspondiente a una lámina de grafeno o a un nanotubo perfectos. 
En realidad, lo anterior no significa realmente que el modelo de campos gauge sea 
incorrecto. Las fases asociadas a —7T3 % 02 y 73 8% 1 tienen su origen en el el hecho 
de que la red es discreta y al empalmado anómalo de la red tras la generación del 
defecto, y son, en principio, imprescindibles en la correcta simulación de la estructura 
electrónica en presencia de estos defectos. En el siguiente capítulo vamos a desarrollar 
un modelo basado en el comportamiento de fermiones de Dirac en espacios curvos, 
en el cual no sólo aparecerá de forma natural la fase asociada a la matriz Y, que 
es la única fase de origen geométrico en el problema, asociada con la conexión de 
espín. Además aparecerán nuevos términos en el hamiltoniano general que originarán 
efectos no triviales, los cuales pueden subsanar las debilidades que presenta el modelo 
de campos gauge. Veremos que en realidad, el modelo de espacios curvos se reduce 
al modelo de campos gauge en el caso de la presencia de un solo defecto, y que es la 
elección de un sistema concreto de coordenadas la que hace que el punto de vista curvo 
se vea enmascarado. Por último, decir que, por supuesto, este modelo de espacios 
curvos admite de forma natural la inclusión de campos gauge externos, por lo que la 
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Figura 3.8: Simulación de uma imagen de microscopía electrónica de efecto tunel 
(STM) en un nanotubo de carbono, extraído de [35]. La diferencia de potencial 
entre la muestra y la punta del microscopio es de 1.5 eV. En los modelos de imagen 
de SIM más simples, la corriente Z que circula entre la muestra y la punta del 
microscopio es directamente proporcional a la densidad de estados local de la muestra: 
Tlr,V) oc p(r,E =V). 


inclusión de los campos gauge asociados a las condiciones de contorno generadas por 
la red se ve salvaguardada. 


Capítulo 4 


Teoría geométrica de los defectos 
topológicos. 





Figura 4.1: Imagen de microscopio de efecto túnel de una membrana suspendida de 
grafeno. En la parte central, señalizada con flechas, se encuentra una monocapa de 
grafeno monocristalina[39]. 


4.1 Introducción. Corrugaciones en grafeno. 
En la introducción del capítulo anterior hemos tratado de dar una visión general 
sobre la descripción de irregularidades y defectos topológicos en membranas cristali- 


nas, suponiendo que tal discusión era aplicable a una lámina suspendida de grafeno. 
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La primera referencia experimental relativa a las corrugaciones en láminas de gra- 
feno suspendidas es la referencia [39]. En este trabajo, por medio de técnicas de 
deposición mecánicas junto con litografía de haz de electrones, los autores fueron ca- 
paces de aislar muestras de grafeno suspendidas, unidas por sus extremos a soportes 
metálicos. 

Las medidas de microscopía de transmisión electrónica han mostrado que en dicha 
configuración, las muestras de grafeno presentan corrugaciones espaciales, tal y como 
se muestra en la figura (4.1). 
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Figura 4.2: Observación experimental de una membrana de grafeno curva[39]. (f): 
Evolución del patrón de difracción con el ángulo de incidencia del haz electrónico. 
Los datos pueden compararse con la simulación de la figura (e). 


Aunque la naturaleza de dichas corrugaciones es aún motivo de debate, trabajos 
recientes [5] muestran que dichas corrugaciones pueden tener su origen en fluctuacio- 
nes térmicas ( recordemos que la teoría general de membranas descrita en el capítulo 
anterior mostraba que las fluctuaciones de la normal a la membrana varían con la 
temperatura como vVT , a temperaturas bajas). 

Otro posible origen de las corrugaciones en las muestras de grafeno es el propio 
desorden morfológico que presenta el sustrato. Experimentalmente, los sustratos so- 
bre los cuales se depositan las muestras de grafeno suelen presentar ona orografía 
compleja, y es de esperar que el grafeno al depositarse, adquiera de forma muy apro- 
ximada la forma del sustrato. Tal efecto se ha comprobado experimentalmente en la 
referencia [40]. En este capítulo vamos a estudiar las propiedades electrónicas de la 
tercera fuente de corrugación en muestras de grafeno: los defectos topológicos. A lo 
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Figura 4.3: (a): Imagen típica de STM de una muestra de grafeno depositada sobre 
un sustrato de SiO). La barra de escala es de 2 nanómetros. (b): Imagen de STM 
con resolución atómica de la misma muestra. La barra de escala es de 2.5 Á. (c) y 
(d): Diferente área de la muestra. La ordenación atómica no es perfecta, indicando 
que existen corrugaciones. 


largo de las páginas siguientes construiremos un modelo basado en las propiedades 
que adquieren los fermiones de Dirac cuando se encuentran en un espacio curvo y 
en presencia de campos gauge externos. Este modelo, que generaliza al modelo de 
campos gauge, nos va a permitir obtener información física para las configuraciones 
de defectos en las cuales dicho modelo se hallaba comprometido, como ocurría con el 
caso de los defectos de tipo Stone-Wales. Además, este modelo de fermiones de Dirac 
ha mostrado su eficacia en la descripción de la estructura de niveles de baja energía 
de la molécula Có0[32, 9]. 

Podemos entender la molécula de Cgy como un fragmento de grafeno en el que se 
han incluido 12 defectos pentagonales, dando como resultado una molécula con forma 
esférica, tal y como se aprecia en la figura (4.1). En las referencias [32, 9] se propuso 
un modelo de espinores en la esfera, en presencia de un campo gauge no abeliano de 
tipo monopolo, que simulaba el efecto de la presencia de los pentágonos en la red 
de grafeno. Las soluciones de la ecuación de Dirac resultante bajo esas condiciones 
se pueden escribir como las autofunciones de un problema de momento angular ge- 
neralizado. Con este formalismo, fue posible deducir los niveles electrónicos y sus 
multiplicidades, permitiendo postular el modelo de fermiones de Dirac en espacios 
curvos y en presencia de campos gauge externos como una potente herramienta que 
puede ser utilizada para estudiar otros tipos de defectos. 
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Figura 4.4: Imagen artística de una molécula de Cép, también conocida como fullereno. 


Por último, en virtud de los resultados experimentales sobre corrugaciones en 
muestras de grafeno, es interesante preguntarse qué tiene que decir el modelo de 
espinores en espacios curvos cuando la superficie que habitan está curvada de manera 
suave, y la red no ha inducido ningún campo gauge externo. Al final de este capítulo 
trataremos de dar respuesta a esta pregunta. 


4.2 Ecuación de Dirac en espacios curvos. 


En la presentación del modelo de espinores en espacios curvos vamos a seguir muy 
de cerca la referencia[41]. 

En general, cuando se trata de trasladar una relación matemática definida en un 
espacio plano a un espacio curvo, simplemente se emplea el principio de covarian- 
cia general, que consiste en sustituir cada magnitud escalar, vectorial o un tensor 
general que se transforma bajo transformaciones definidas en ese espacio plano por 
su homólogo, que se transforma bajo transformaciones generales en el espacio curvo. 
Esta simple sustitución no funciona para objetos que se comportan como espinores 
bajo transformaciones de coordenadas. La razón matemática es que existen represen- 
taciones tensoriales del grupo de transformaciones generales, pero no existen repre- 
sentaciones espinoriales de dicho grupo. Esta consideración hace que haya que buscar 
una formulación alternativa para introducir espinores en un espacio curvo. 

El formalismo adecuado para el estudio de espinores en espacios curvos es el 
formalismo de las tétradas. El formalismo de tétradas se basa en el Principio de 
Equivalencia, que afirma que en cada punto X de una superficie curva, es posible 
definir un sistema de coordenadas £x que es localmente plano en dicho punto. En 
general, dicho sistema de coordenadas no será el mismo en todo punto de la superficie 
a menos que la superficie sea realmente plana. En general, podemos escribir el tensor 
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métrico en dicho sistema de coordenadas como 





_ a,b 
9uw = € €, Mab» (4.1) 
a 0% : e 
donde ef, = ¿2 se transforman como un conjunto de vectores covariantes bajo un 
cambio de coordenadas general: 
Ox” 
a la a 
> la = tr (4.2) 


Entendiendo este conjunto de vectores como una base en el espacio vectorial (co- 
tangente) definido en cada punto, podemos proyectar cualquier vector contravariante 
A*(x) en el punto x sobre dicho conjunto de vectores: 


+A" =e2AP, (4.3) 


Podemos generalizar la relación anterior a vectores covariantes y a tensores en general, 
definiendo los vectores e. Nótese que los índices latinos bajan y suben empleando la 
métrica plana Nay y su inversa respectivamente, y los índices griegos hacen lo mismo 
pero con la métrica curva 9y,. Es inmediato comprobar que, con estas reglas, el tensor 
métrico en estas coordenadas es, como era de esperar, *9ap = Nab- 

Las cantidades *A“ constituyen un vector contravariante en el espacio plano defi- 
nido en el punto X, por lo que podemos tomar la alternativa de olvidar las cantidades 
tensoriales definidas en un espacio curvo general, y trabajar con el sistema de coorde- 
nadas localmente plano en cada punto que determinan ef, y las cantidades tensoriales 
definidas en ese espacio plano. 

Una vez que estamos en el sistema de coordenadas plano, cualquier otro sistema 
de coordenadas plano que se pueda definir será igual de válido que el inicial, por lo 
que podremos transformar, en el punto escogido x a través de una transformación de 
coordenadas plana A(x), las cantidades tensoriales planas xB%-",,,...: 


BAN BA (4.4) 


Esta regla de transformación se puede generalizar para objetos que se transformen 
como espinores bajo cambios de coordenadas planos: 


donde [D(A)] nm es una representación espinorial de la transformación A(x). 

Para escribir el hamiltoniano de Dirac en una superficie curva, debemos saber 
también cómo se transforma la derivada de un espinor. La derivada es un objeto que 
se comporta como un vector bajo cambios generales de coordenadas: 


0 0%" 0 


Oxu "Damm dar 





(4.6) 
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Como en los ejemplos anteriores, podemos proyectar dicho vector sobre el conjunto de 
vectores (e) y, una vez escrito en la base plana, aplicarlo sobre una magnitud espino- 
rial. La derivada se transformará entonces bajo una transformación de coordenadas 
plana A(x) como 


AaJleh > [D(A(2)) +0(2)] = (4.7) 


= loli [IA 0 + (¿GDA +00) 


Como podemos ver, la derivada de un espinor no se transforma como un vector 
bajo cambios de coordenadas planas. Este problema se resuelve definiendo un vector 
derivada covariante para espinores que tenga la forma 


0) 
Da = e fa + r, A (4.8) 
y que se transforme de forma usual bajo transformaciones planas tanto en el índice 
vectorial plano como en el índice espinorial: 


Da * Ya — MD(NamD; * Um. (4.9) 


Las matrices I', han de estar definidas de tal forma que su ley de transformación bajo 
cambios de coordenadas planas sea: 


0 


PDD) [72 


Dn) D(A). (4.10) 
El segundo término en (4.10) va a cancelar el segundo término en (4.8) para dar lugar 
a la ley de transformación deseada (4.9). 

Para obtener la estructura de las matrices I',(x) es suficiente considerar la manera 
en la que se comportan dichas matrices y los vectores Le? bajo transformaciones en 
el espacio plano. El resultado es que la ley de transformación (4.10) es satisfecha por 
las matrices 


1 a v 
Ls e ATA (4.11) 


donde e? es la derivada covariante de los vectores ef: 
bip db 
Po, v vo p 


Por lo tanto, empleando este formalismo, la forma covariante de la acción en un espa- 
cio curvo se obtiene reemplazando cada magnitud B”--,,,, por su correspondiente en 
el sistema de coordenadas planas locales xB"*--,,,,..., y las derivadas por las derivadas 
correspondientes dependiendo de si la magnitud es un vector (por ejemplo (4.12)) o 
una magnitud espinorial (4.8). 
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En concreto, el hamiltoniano para un espinor de Dirac sin masa en un espacio 
curvo general se escribe de la siguiente manera: 


HS Sa RD, +T,)4(x). (4.13) 


En la expresión (4.13), aparte de la sustitución de la derivada covariante para el 
espinor 4 (4.8), las matrices y(x) dependientes de x están definidas a través de (4.3) 
como 
yx) =*y"(x) = y eq (1). (4.14) 

El hamiltoniano (4.13), de esta manera, es invariante bajo transformaciones de coorde- 
nadas generales, independientemente de la elección concreta de sistema de referencia 
plano £% escogido en cada punto. 

En las siguientes secciones veremos de qué manera el formalismo de espacios curvos 
completa el modelo de campos gauge descrito en el capítulo anterior, y veremos que 
información adicional nos proporciona dicho formalismo. 


4.3 Descripción de un número arbitrario de discli- 
naciones. 


Comenzaremos esta sección conectando el lenguaje de espacios curvos con el mo- 
delo para una disclinación aislada, descrito en el capítulo anterior. 

Para tal fin, el primer paso es encontrar el tensor métrico que define una superficie 
cónica, que es la superficie que resulta de eliminar un sector angular de superficie plana 
(ver, por ejemplo, la figura (3.4). El tensor métrico en la superficie plana se escribe, 
en coordenadas polares planas como: 


ds? = dr? + r?d6?, (4.15) 


donde las variables r y 0 están definidas en (0,00) y (0,27), respectivamente. En 
el caso de una superficie cónica, la variable angular O tendrá un dominio menor, 
concretamente, 0 € (0,27 — 5), siendo 8 el déficit angular que posee el cono con 
respecto a la superficie plana. Con el simple cambio de variables (r,0) = (r”,0”) = 
(r, (1— 4)0) podemos transformar el tensor métrico definido a partir de (4.15) a 


2 
ds? = dr? + ( = Z) r?de?, (4.16) 
T 


por lo que las componentes no nulas del tensor métrico son 


B 2 
9rr = 1, 906 — ( A Z) r?. (4.17) 
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A partir de ahora, por sencillez, definiremos a = (1 — 2). Utilizando el formalismo 


descrito en el apartado anterior, podemos calcular el conjunto de vectores Le? o 
díadas y la conexión de espín. En coordenadas polares, las matrices y* y la conexión 
de espín tienen la forma: 


. 1 
=> Y =—Y, (4.18) 
QT 


T,=0, To =-an. (4.19) 


El hamiltoniano en coordenadas polares (4.13) se escribe entonces como 
As ES ivp fa ry a (ad, E 2) eN 00] Y. (4.20) 
F r 


Retrocedamos brevemente al capítulo anterior, concretamente a las ecuaciones (3.14) 
y (3.17). De esas expresiones podemos inferir la forma funcional que tendrá un campo 
gauge general que procede de una condición de contorno, ya que la única manera para 
que se verifique $ Adr = P es que A(r) tenga estructura de vórtice: 


(4.21) 


En coordenadas polares, el campo A(r) sólo tendrá componente angular: A, =0 y 
Ay = —. Puesto que en la ecuación (3.19) el campo gauge asociado a 2 entra en forma 
de producto escalar, iy" A;, siendo el producto escalar un invariante bajo cambio 
de coordenadas, en coordenadas polares tendrá la forma ¡y?L.Ajy. Como Y = y? 
yq La?, concluimos inmediatamente que, en coordenadas polares, el término 


asociado al campo gauge en el hamiltoniano (3.19) es Ap = pl. 
Podemos escribir el hamiltoniano en coordenadas polares planas como 


- 1 
Hue + iop | dro y (o, == 2) + 17) Y (4.22) 


y compararlo directamente con (4.20). La primera apreciación es que en ambos for- 
malismos, como ya hemos mencionado anteriormente, aparece el mismo término co- 
rrespondiente al campo gauge. Sin embargo, existe una diferencia entre ambos ha- 
miltonianos, correspondiente al término de la derivada radial: mientras que en (4.20) 
el término es ay*0,, en (4.22) dicho término es simplemente y*0,.. 

Podemos interpretar esta diferencia en términos de la aparición de una velocidad 
de Fermi modificada. Consideremos de nuevo el hamiltoniano de Dirac en una su- 
perficie plana, y calculemos el operador velocidad en primera cuantización. De las 
ecuaciones de Heisenberg, tenemos que 


v” =-—¿|Hp,x*] = upy". (4.23) 
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En coordenadas polares, y en el espacio plano, el operador velocidad tiene como 


componentes 


1 
VU, = Up Vo = UY, (4.24) 


sin embargo, en las coordenadas definidas en (4.16), con el hamiltoniano definido en 
(4.20) dicho operador velocidad tiene la expresión 


Ccuru I 


ver = avpy! DE UY, (4.25) 


comparando (4.24) y (4.25) llegamos a la conclusión anunciada anteriormente: en 
el caso de una disclinación aislada, hemos convertido el hamiltoniano de Dirac curvo 
(4.13) en el hamiltoniano plano con un campo gauge (4.19) y una velocidad modificada 
(4.25), ausente en el modelo de campos gauge (4.22). La aparición de una velocidad 
modificada por la existencia de defectos (y en general, por el hecho de tener espinores 
inmersos en un espacio curvo) hace que el modelo de espacios curvos sea más rico 
que el modelo de campos gauge, pudiendo obtener variaciones en las propiedades 
electrónicas respecto a la situación ideal en configuraciones de defectos en las cuales 
el modelo de campos gauge no ofrecería nada, como es el caso, por ejemplo, de 
los defectos Stone-Wales. En (4.22) y (4.20) no se han incluido los campos gauge 
correspondientes a las condiciones de contorno cuyo origen es la red. La forma para 
incluirlos es semejante al campo (4.21), puesto que en todos los casos el campo gauge 
es singular, solamente cambia la matriz que acompaña a dicho campo. La ecuación 
(4.13) para espinores en una superficie curva que estén acoplados a un campo gauge 
general A,,, tiene la forma 


Hito —= 5 aryivajior(a) (O, en ES as ¿Qa Ay a) (2). (4.26) 


Nuestro objetivo es entonces aplicar el formalismo de espacios curvos cuando con- 
sideramos la situación de una lámina de grafeno que contiene un número arbitrario de 
defectos localizados en posiciones conocidas. Para tal fin, procederemos de manera 
similar al caso de una disclinación aislada, encontrando una métrica que describa 
apropiadamente el espacio curvo generado por un conjunto de disclinaciones, y cal- 
culando posteriormente los vectores ef y la conexión de espín. 

Podemos trabajar en este caso por analogía. Las disclinaciones en grafeno son el 
análogo bidimensional a las cuerdas cósmicas en 3 +1 dimensiones|42, 25, 43]. 

Podemos encontrar la métrica que describe un conjunto arbitrario de defectos a 
partir de la métrica para un defecto aislado (4.16) empleando los siguientes cambios 
de variables: 

P=iane 8=0 (4.27) 


=P co, y =r'sin(0”. (4.28) 
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La métrica (4.16) se escribe entonces, en coordenadas conformes (x, y), como 
ds? = e2A(r) (da? + dy?) , (4.29) 
quedando definido el factor conforme A(zx, y) de la siguiente manera: 
Me y) = og: (4.30) 


quedando el parámetro A definido como A= 1 -— «a = Z. Para el caso de un número 
N de defectos situados en las posiciones r;, y con una amplitud angular arbitraria Ay, 
generalizamos el factor conforme (4.30) a 


Ma 2 A¡log(|r' — r/)). (4.31) 


Podemos calcular, a modo de comprobación de la generalización (4.31), la cur- 
vatura intrínseca correspondiente a la métrica (4.29). En coordenadas conformes|44] 
la curvatura intrínseca para una superficie bidimensional descrita en coordenadas 
conformes satisface la siguiente relación 


R=e2ArWy?p, (4.32) 


que con (4.31) da R = 450(r — r;). En coordenadas polares resulta ser la expresión 
conocida 


R= y da (r —r;), (4.33) 


es decir, la curvatura generada por los defectos sobre la superficie bidimensional 
está concentrada en los sitios donde se encuentran los defectos, siendo el resto de la 
superficie plana. 

Escogeremos el sistema localmente plano definido a partir de los vectores ez, que 
tengan las siguientes componentes el = e? =e714 y e) = e? =0. En dicho sistema de 
coordenadas localmente plano, la conexión de espín es 


1 
Fis 31104, (4.34) 


1 
> = ¿11 02h. 


4.4 Densidad de estados. 


A continuación analizaremos las consecuencias físicas que conlleva la existencia de 
defectos topológicos en una lámina de grafeno, en concreto, analizaremos los cambios 
en la densidad local de estados. En los ejemplos que discutiremos a continuación 
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hemos optado por mantenernos en el nivel de complejidad más bajo, de tal forma 
que estudiaremos configuraciones de defectos para las cuales las contribuciones de los 
campos gauge asociados a las matrices % y Rp son nulas, tal como discutimos en el 
capítulo anterior. De este modo, podremos centrarnos en entender de qué manera la 
curvatura de la lámina de grafeno modifica la estructura electrónica dentro de este 
modelo geométrico. 

Como sabemos del capítulo 1, la densidad de estados, definida como el número 
de estados electrónicos por unidad de área y de energía, se puede calcular a partir de 
la función de Green a dos puntos, siendo la densidad de estados la parte imaginaria 
de dicha función de Green. En una lámina de grafeno perfecto, debido a que existe 
simetría traslacional, dicha densidad de estados es independiente de la posición y es 
proporcional al valor absoluto de la energía. En presencia de desorden, la simetría 
bajo translaciones se pierde y la función de Green a dos puntos deja de depender de 
la combinación r — r' y por lo tanto la densidad de estados pasará a depender de la 
posición en la lámina: 

p(r, E) = lim, TrImG(E,r,r My. (4.35) 
En lo que sigue, vamos a calcular la función de Green a dos puntos resolviendo la 
ecuación de Green en espacios curvos: 


VEG(E,r,r') —ivpy(x) (9, +T,) G(E, r,r”) = OS —r'). (4.36) 


La ecuación (4.36) es en principio, demasiado complicada para resolverla de forma 
general, por lo que optaremos por resolverla de forma perturbativa en el parámetro 
A obteniendo correcciones perturbativas a la densidad de estados local. 

Desarrollando en serie de potencias de A en (4.36), obtenemos a primer orden en 
A la siguiente ecuación: 


VEG(E,r,r) — iv 0,G(E, r,r') —- V(E,r)G(E,r,r”) =9% (r — r), (4.37) 


que interpretamos como la ecuación de Green para la función de Green en el espacio 
plano en presencia del potencial externo V(£, r), 


V(E,r) = —EA(9) + ivpA(m)y?0; + E (0,A(r)). (4.38) 


En la expresión (4.38) reconocemos términos similares de velocidad dependiente de 
la posición generada por el desorden 


ivpA(r)y?0, 
que aparece en (4.25) y el término correspondiente al campo gauge 


1Up 


y? (01) 
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, análogo a (4.21). También obtenemos un término extra, dependiente de la energía, 
que no aparece en el análisis que hicimos sobre el hamiltoniano. Este término es 
consecuencia de la estrategia que estamos utilizando para calcular la densidad de 
estados a partir de la resolución de la ecuación de Green (4.36). Una vez que hemos 
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Figura 4.5: Representación pictórica de la ecuación integral (4.39). La línea vertical 
continua representa a G(E,r — r') y la estrella representa el potencial V(£, r;). 


transformado el problema de resolver la ecuación de Green en un espacio curvo al 
problema de resolver la ecuación de Green en un espacio plano con un potencial de 
desorden V(E, r), podemos emplear la metodología estándar diseñada para resolver 
estos problemas[45]. 

La ecuación diferencial (4.37) puede convertirse en la siguiente ecuación integral 


G(E,r,r') =Go(E,r-1)+ / BroGo(E,r—ra)V(E,ro)JG(E, ro, 10), (4.39) 


que resolveremos de forma iterativa en serie de potencias de A, parámetro que se 
encuentra incluido en la función A(x, y). Por consistencia, nos quedamos en el primer 
orden en A, ya que en la deducción de la ecuación (4.37) nos quedamos con el primer 
orden en A. El cálculo de la función de Green a mayor orden en A requeriría calcular 
V(E,r) a orden O(A?), así como el término de orden O(A?) para el potencial (4.38). 
A orden O(A), la función de Green (4.39) se escribe como 


G(E,r, 1) = Go E,r-r)+ a drsG(E,r —ro)V(E,ro)JGo(E,ra Y). (4.40) 


Puesto que se trata de utilizar la función de Green para calcular las correcciones a 
la densidad de estados Ap(E,r) = p(E,r) — p"(E,r), nos centraremos en el segundo 
término de (4.40), del cual extraeremos Ap(E, r). 
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Empleando la expresión de la función de Green libre: 


2% / EY—opykiX 
Go(E,r 3% r”) > d Y UFY k; ) e k(r=r) (4.41) 








4712 | E2 — vhk? + ie 
y si 
= np 1pr> 
Ara) = [e*A(p), (4.42) 
con 
Ed Ajerti 
A(p) = Y, 7 (4.43) 
J 


en el límite r' > r la traza del segundo término de la ecuación (4.40) queda 
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Figura 4.6: Izquierda: Imagen de la corrección a la densidad local de estados inducida 
por un par heptágono-pentágono en una porción amplia de la muestra de grafeno, a 
una energía normalizada w = 0.7. El par de defectos se encuentra localizado en el 
centro. La imagen muestra que la corrección a la densidad de estados es mayor en 
puntos próximos a las posiciones de los defectos. Derecha: Detalle de la densidad de 
estados con el par localizado fuera de la imagen. 


TrG(B,r,r) = 5 57z el pe A(PT(E, p), (4.44) 


donde T(E, p) es 
8E3 — 2Eg? 


T(E, sl a a 

ERE CAN 1) 

La corrección a la densidad de estados Ap(E,r) se obtiene, como hemos dicho, to- 
mando la parte imaginaria de (4.44): 


DEAL Jo (Er; 
AE, Y) =>3 al z / a apo —F(E,k), (4.46) 





(4.45) 
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donde Jo(Er;) es la función de Bessel de primera especie de orden 0, y la función 
F(E, k) viene dada por las siguientes expresiones 





2-2 
- arg tam il 


Y) 4H? < k2 
F(E,k)= VTA (4.47) 


t _ 
OA EVITE AE? > p?2 
V4E2-—k2 p) sz 








A continuación utilizaremos el resultado (4.46) para discutir algunos casos de 
interés físico. Comenzaremos por la situación en la que tenemos una única pareja 
de defectos pentágono-heptágono situados en posiciones (—1,0) y (1,0), en notación 
quiral. En la figura (4.6) se muestra la corrección a la densidad de estados para esta 
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Figura 4.7: Corrección a la densidad de estados en diferentes puntos de la muestra, 
indicados en la leyenda, como función de la energía. El par heptágono-pentágono se 
encuentra localizado en las mismas posiciones que en la figura (4.4). 


disposición de defectos a una energía (normalizada a una frecuencia de corte (Q) de 
w= 0.7. El color verde corresponde a las zonas donde la corrección es casi nula o nula, 
el color rojo corresponde a zonas donde aparece acumulación de densidad de estados, 
siendo el azul correspondiente a zonas donde la densidad de estados ha disminuido. 
En la figura (4.6, izquierda) la imagen está dominada por la densidad de estados en 
las regiones que rodean ambos defectos. En torno a la posiciones de los defectos, la 
corrección a la densidad de estados (4.46) es divergente, y la validez de la aproximación 
(4.40) se pierde. Sin embargo, lejos de las posiciones de los defectos, donde las 
aproximaciones realizadas tienen validez, observamos una modulación de la densidad 
de estados que presenta un marcado carácter bipolar, hecho que está de acuerdo 
con simulaciones numéricas|46). La figura (4.7) muestra la densidad de estados para 
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Figura 4.8: Densidad de estados local para dos pares de defectos heptágono-pentágono 
para cuantro energías diferentes. 


la misma disposición de defectos que la figura (4.6) como función de la energía, en 
distintos puntos de la lámina. En esta figura se ve de forma clara que la densidad de 
estados presenta oscilaciones a medida que la variamos la energía, y dichas oscilaciones 
dependen fuertemente de la posición sobre la lámina de grafeno donde se observen. 
En una muestra real, los defectos se encuentran repartidos de forma aleatoria y, 
obviamente, en mayor número, así que podemos estudiar configuraciones de defectos 
que, cumpliendo las ligaduras impuestas para la anulación de los campos gauge extras, 
puedan representar esta situación más realista. En la figura (4.8) mostramos, para 
diferentes energías, la corrección a la densidad de estados local para dos defectos 
compuestos de un heptágono y un pentágono. Á energías bajas, la corrección es casi 
cero, salvo en las regiones que rodean los defectos, en donde, como hemos dicho, la 
aproximación adoptada carece de validez. A medida que aumentamos la energía, la 
densidad de estados se va viendo modificada en forma de oscilaciones inhomogéneas, 
debidas principalmente a la orientación relativa de los defectos (nótese que se sigue 
apreciando el carácter dipolar en las proximidades de los pares). También es relevante 
mencionar que no parece que la longitud de onda de las oscilaciones escale de manera 
inversa con la energía, como ocurre con un dipolo en electromagnetismo clásico. 

Los defectos de tipo Stone-Wales constituyen otra disposición de defectos de in- 
terés para estudiar. Como ya sabemos, estos defectos son los que tienen mayor pro- 
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Figura 4.9: Corrección a la densidad de estados local para un defecto de tipo Stone- 
Wales. 








QU ON >= oO >. u.ée 


ES 


babilidad de ocurrir en la realidad, ya de forma natural o en muestras artificialmente 
dañadas. 

En la figura (4.9) se muestra la densidad de estados a una energía normalizada 
constante w = 0.7 en una muestra amplia de grafeno que contiene en su centro un 
defecto Stone-Wales. Como en los casos anteriores, las modificaciones en la densidad 
de estados se encuentran fuertemente atenuadas, debido a la presencia de las correc- 
ciones no válidas en los entornos muy próximos al defecto. En la figura (4.9, derecha), 
se muestra una sección de lámina de grafeno con el defecto localizado fuera del plano. 
Como antes, se observan claramente ondulaciones ordenadas en la densidad de es- 
tados, debidas a la simetría del defecto, y pueden compararse, al menos de forma 
cualitativa, con modificaciones a la densidad de estados encontradas en simulaciones 
de nanotubos que contienen dichos defectos[35). 


4.5 Efecto de la curvatura suave en la estructura 
electrónica. 


En las secciones anteriores hemos empleado el formalismo de espinores en espacios 
curvos para describir el comportamiento electrónico de baja energía en una lámina 
de grafeno con curvatura inducida únicamente por defectos topológicos. 

En la introducción del presente capítulo hemos visto que gran la mayoría de mues- 
tras de grafeno depositadas sobre un sustrato[40] o suspendidas[5] presentan corruga- 
ciones espaciales, visibles a través de microscopía de por transmisión de electrones[47] 
y por espectroscopía de efecto túnel[48]. 

Aunque la existencia de defectos topológicos es una fuente de corrugación en 
grafeno, no podemos decir que sea la única: fluctuaciones térmicas y corrugaciones 
en el sustrato inducirán curvatura en una lámina de grafeno. Por otro lado, sabemos 
que la curvatura inducida por los defectos es singular y está localizada donde se sitúa 
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el defecto. La pregunta que nos hacemos ahora es qué efectos podemos esperar de 
una muestra de grafeno donde, sin especificar el motivo, se presentan regiones con 
curvatura suave y más o menos extensa. 

Para responder esa pregunta contamos con las herramientas matemáticas que 
hemos descrito en las secciones anteriores, por lo que seguiremos el mismo programa 
que propusimos para el estudio de los efectos en la estructura electrónica en grafeno 
curvado por defectos topológicos. 

El primer paso es describir la superficie de grafeno curvado a través de una métrica. 
Consideraremos que la superficie, inmersa en el espacio tridimensional, se puede des- 
cribir por medio de la función z(r) que representa la altura de la superficie con respecto 
a la situación plana z = 0. Por sencillez, consideraremos regiones curvas con simetría 
bajo rotaciones, por lo que z será función únicamente de la coordenada radial r: 
z(r) = z(r). La métrica de la superficie vendrá descrita entonces por el elemento de 
línea 

ds? = dr? + r?d6? + dz? = (1+af(r)) dr? + r?d8*?, (4.48) 


donde hemos escrito (a es un parámetro libre que emplearemos para cuantificar cuánto 
nos alejamos de la situación plana) 


de? = pra dr? = af(r)dr?. (4.49) 


Los únicos requisitos sobre la función f son que la superficie sea asintóticamente 


100 





Figura 4.10: Forma de la corrugación descrita por la función z(r) para simular la 
superficie curva suave de grafeno. 
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plana, es decir, f(r) — 0 en el límite r —> oo, y que la superficie z(r) sea suave en el 
origen r = 0, condición que se traduce en que la función f se anule suficientemente 
rápido en el origen. El ejemplo más sencillo de dicha situación es el de una corrugación 
con forma de gaussiana, que utilizaremos para mostrar de qué manera funciona el 
modelo: 


2 
z= Áexp (5) a (4.50) 


donde A es la altura máxima de la gaussiana, y b es la anchura de la gaussiana. Con 
esta elección para la función z(r) la función f(r) y el parámetro a: quedan definidos 


— e (2) exp (37) (4.51) 


y a = (4/b)?. Las corrugaciones que se observan experimentalmente nos proporcio- 
nan valores típicos para el cociente entre la anchura y la altura de la gaussiana se 
pueden encontrar en las referencias [39] y [49], variando en un rango de entre 0.1 y 
0.3. Dicho rango se traduce para el parámetro a en 0.01 y 0.1, aproximadamente, 
permitiéndonos estudiar el efecto de la curvatura en la densidad de estados en el 
régimen perturbativo del parámetro «a. 

Los vectores ez, calculados para la métrica definida en (4.48) se pueden escribir de 
forma matricial como 


pe Laf?coso —rsin0 


(0-7 
n= (1+af] sing rcosQ la 


Los vectores descritos en (4.52) no se reducen a los vectores ez, planos cuando a: = 0 
que corresponderían a la situación plana en (4.18). Esto se debe a que la condición 
(4.1) no define de forma unívoca los vectores ef, y la elección de un conjunto u otro 
no tiene ninguna relevancia a la hora de calcular magnitudes físicas. Este hecho va a 
quedar patente una vez conocidas las componentes de la conexión de espín: 


1-(1-ap 2 169 
9 Y 


La elección de díadas (4.53) no se reduce al caso plano descrito en la sección anterior, 
cuando estudiamos el efecto de una disclinación aislada en coordenadas polares en el 
límite plano 6 = 0. Esto se debe simplemente a la libertad definición de las díadas 
que nos permite la ecuación (4.1). La ley de transformación de (4.52) en el límite 
a =0 al sistema de vectores ey, descrito en la sección anterior en el límite 6 =0 


O 
Á=(. »), (4.54) 


es una simple rotación de ángulo 9 en torno al eje perpendicular a la superficie: 


ma) = cosó sinó ly 


(4.52) 





T,=0, Ty = (4.53) 





—sinó cosó (4.55) 
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La representación espinorial (4.5) de la transformación (4.55) es la conocida matriz 


D(A*,) = exp (597) (4.56) 


Introduciendo esta D(A) en la expresión (4.10) obtenemos la conexión de espín para 
el sistema de coordenadas decrito por (4.54) en el límite plano: 
8 0 11 
Py =9%/0,0(1) =>A9= y a 
como cabía esperar. 

Para calcular la densidad de estados, seguiremos la estrategia definida en la sección 
anterior para el caso de múltiples defectos. Una vez calculada la conexión de espín 
y las díadas, resolveremos perturbativamente la ecuación de Green (4.36) para este 
caso, y de ella extraeremos las correcciones a la densidad de estados local. Aunque la 
estrategia a seguir sea la misma, vamos a proceder de un modo ligeramente diferente 
al caso anterior, para aprovechar al máximo las ventajas de la simetría rotacional que 
presenta el problema. 

Puesto que vamos a trabajar a orden O(a), desarrollaremos el término y que 
multiplica a 4% (r — 1”) en el lado derecho de (4.36): g71? = 1(1— “£). En cuanto 
al lado izquierdo de (4.36), el desarrollo en serie a orden O(«a) genera el término de 
potencial 


-1/2 


2 2r 
en donde la función I'(9) quedará definida más adelante. En el potencial (4.57) junto 
con el término libre que aparecerá en la ecuación de Green, podemos identificar los 
términos de velocidad de Fermi modificada en la dirección radial, a orden O(a): 


Vr, 0) = PESO pg) C -0,), (4.57) 


o =0pP(O/1— 5aJ(r)), 


y el potencial magnético 





El campo magnético asociado a Ag, a orden más bajo en el parámetro a es 
B.=af(r) (4.58) 
¿= —af'(r). , 
dr 


La magnitud de este campo magnético resulta ser del orden de 0.5 a 2 Teslas, compa- 
tible con las estimaciones dadas en la referencia [49], siendo una posible explicación a 
la ausencia de localización débil descrita allí. Es interesante notar que la velocidad de 
Fermi modificada por la curvatura es siempre menor que el valor vp, independiente- 
mente del tipo de superficie considerada. Inspiradas por este trabajo se han realizado 
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recientemente medidas de microscopía túnel sobre muestras de grafeno suspendido 
que sugieren la existencia de una variación local de la velocidad de Fermi. El análisis 
de estos resultados experimentales en el marco de la teoría de espinores en grafeno 
curvo será objeto de un estudio futuro. 

A continuación expondremos los resultados obtenidos para las correcciones a la 
densidad de estados. El primer término en el desarrollo de G(E, r, r”) equivalente a 
(4.40), teniendo en cuenta los dos términos originados por g”*/? y la representación 
de la función delta de Dirac en coordenadas polares, es: 


1 
G(E,r,r') = Go(£, r=r)-0f(r)Go(E, r—r)+f ProGo(E, r—r2)V (r2)Go(E, rar”). 
(4.59) 
Aunque en el caso de múltiples defectos topológicos, trabajamos con la representación 


o 


g 





Figura 4.11: Corrección a la densidad de estados local bajo la corrugación en una 
lámina de grafeno. Color oscuro significa una contribución negativa a la densidad 
de estados, mientras que la escala de colores clara indicará contribuciones positivas 
respecto de la densidad de estados libre para el grafeno a la energía dada. 


de momentos del propagador, aquí nos mantendremos en el espacio real, por lo que 
tendremos que definir el propagador de Dirac (4.41) en coordenadas polares. La forma 
más sencilla es aplicar la bien conocida relación existente entre el propagador de Dirac 
y el propagador de Klein-Gordon 


Cubiras => DO Curia Gdor (4.60) 
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a través del operador de Dirac en coordenadas polares: 
DP (rio, q 1(9)20,) (4.61) 
donde 
T(0) = y cos0 + y sinó, (4.62) 
T"(9) = —y! sin 0 + y cos60. 


El propagador de Klein-Gordon en el espacio de posiciones tiene la forma 
1 
Goktem=Gordon (Le, Ls r”) + >= HolElr NN r)), (4.63) 


siendo Ho la función de Hankel de orden cero(50]. La aplicación del operador (4.61) 
a (4.63) genera la siguiente expresión para el propagador de Dirac en términos de Hoy 
y la función de Hankel de orden uno: 


=1B E 


GopiradE,r — Y”) = — VHlE [lr — r”|) — a. e lr —r)) (rT (9) — r'T (05). 


(4.64) 
Usando las expresiones (4.64), (4.59) y (4.35), llegamos a la expresión siguiente para 
la densidad de estados local: 








; 1 af(r) a 
p(E,r') do aria 5 e J aras 
, r —r'cos(0 — 0”) 
e a 
2 () + rf Y (EAN (EAN), (4.65) 


en donde hemos definido Ar = |r —7”| y Jo 1, Yo 1 son funciones de Bessel de primera 
y segunda especie, respectivamente. Exceptuando el cálculo explícito de los factores 
geométricos, en los cálculos anteriores, no se ha hecho mención explícita a la forma 
concreta de la corrugación, por lo que, aunque en las figuras se haya escogido una 
forma de gaussiana, la discusión que sigue será general para cualquier tipo de función 
z(r) axialmente simétrica. 

De la ecuación (4.65) se aprecia que, debido a la presencia de funciones de Bessel, 
la corrección a la densidad de estados muestra oscilaciones cuya frecuencia crece al 
aumentar la energía y cuya amplitud decrecerá con la distancia como r”!. En la figura 
(4.11) se muestra la corrección a la densidad de estados local para una corrugación 
correspondiente a los parámetros b= 504 y una energía E = 0.1eV. El máximo valor 
de la corrección con respecto al valor de la densidad de estados libre a esa energía 
es del orden de un uno por ciento, para una corrugación correspondiente al valor 
a “0.01. 
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Figura 4.12: Comparación entre la corrección a la densidad de estados (linea de 
puntos) en unidades arbitrarias, la forma de la corrugación (color rojo) y la curvatura 
de dicha corrugación (color amarillo). 


En la figura (4.12) comparamos la corrección a la densidad de estados con la mor- 
fología de la corrugación (en este caso sí definida como una gaussiana) y la curvatura 
de la superficie descrita por dicha gaussiana. A bajas energías, se puede ver que el va- 
lor mínimo de la corrección está correlacionado con el el primer cero de la curvatura, 
que, de nuevo para el caso particular de la gaussiana, corresponde a v/2b. 


En el caso general de una corrugación de perfil arbitrario, se puede demostrar que 
el potencial gauge efectivo no contribuye a la corrección a la densidad de estados a 
primer orden en teoría de perturbaciones. En este caso, como el caso de los defectos 
topológicos, el campo inducido será proporcional a la matriz y?, lo que significa que 
de manera efectiva, cada uno de los dos conos en el la relación de dispersión notará un 
campo magnético de signo contrario, manteniéndose el sistema completo invariante 
bajo inversión espacial. En equilibrio, la corriente total será cero, lo que significará 
que la carga no se acopla al campo magnético, al menos en primer orden de teoría 
de perturbaciones[51]. Esta consideración explica el hecho de la anulación de la 
contribución a la densidad de estados asociada al campo gauge. 


Por último, comentaremos que el valor máximo de la corrección a la densidad de 
estados se concentra en torno al tamaño típico de la corrugación, por lo que en un 
experimento de STM, las modificaciones más altas a la densidad de estados deberán 
estar en correspondencia con la morfología de la superficie. 
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4.6 Conclusiones. 


En el presente capítulo, hemos planteado un formalismo alternativo al modelo de 
fermiones de Dirac en presencia de campos gauge, estudiado en el capítulo anterior. 
Dicho formalismo, consistente en el estudio de dichos fermiones en un espacio curvo 
junto con los campos gauge mencionados generaliza y complementa el modelo de 
campos gauge, pudiéndose analizar disposiciones de defectos para las cuales el mero 
modelo de campos es sospechoso. Hemos calculado perturbativamente las correcciones 
a la densidad de estados local para varias disposiciones de defectos, algunos de ellos 
de relevancia física, como el defecto Stone-Wales, para el cual el modelo de campos 
gauge no da corrección alguna. 

Las corrugaciones en grafeno pueden tener más orígenes aparte de los defectos 
topológicos. En la última parte de este capítulo, hemos aplicado el modelo de espacios 
curvos a la situación de una lámina de grafeno con curvatura suave. Hemos calculado 
las correcciones a la densidad de estados válidas para corrugaciones con simetría axial, 
aunque hayamos tomado como referencia una corrugación con forma gaussiana. Junto 
con oscilaciones características, hemos encontrado que el campo gauge asociado a la 
curvatura no contribuye a estas correcciones de la densidad de estados, al menos a 
primer orden en teoría de perturbaciones, como era de esperar. 


Capítulo 5 


Desorden topológico y transporte 
en grafeno. 





5.1 Introducción. Situación experimental. 


En muestras experimentales de tamaños típicos del orden de unas cuantas micras, 
el desorden, independientemente de su origen, puede jugar un papel relevante en 
algunas de las propiedades físicas del sistema. En muestras de ese tamaño, la variedad 
de formas de distribuirse dicho desorden sobre la muestra es enorme. Esto hace que 
cuando se mida sobre diferentes muestras, encontremos magnitudes que varíen de 
muestra a muestra de manera aleatoria; magnitudes que, aún variando entre muestras, 
presenten un patrón que puede ser estudiado en forma de ley física; y por último 
magnitudes que tienen un caracter universal, independiente de cual muestra se esté 
analizando. Para estas dos últimas categorías hay que establecer un modelo para el 
desorden en el que la distribución concreta de defectos sobre la muestra sea irrelevante. 
A menudo el desorden puede describirse por medio de un potencial V (r,r,) que 
depende explícitamente de las posiciones r, de los defectos. El formalismo a utilizar 
debe permitirnos calcular observables independientes de r,, pero dependientes de unos 
pocos parámetros que caractericen el tipo de desorden, y eventualmente, la densidad 
de defectos presente. 

Desde las primeras observaciones experimentales sobre conductividad realizadas 
sobre muestras de grafeno[52, 53], una de las propiedades más sorprendentes y que 
más atención ha recibido es la obtención de un valor mínimo finito, del orden del 
cuanto de conductancia — 4e2/h. La observación más robusta en situación de medio 
llenado parece ser independiente de la temperatura en un amplio rango de tempera- 
tura. Otro hecho a destacar es la dependencia de la conductividad con el potencial 
de puerta aplicado, o lo que es lo mismo, con la densidad de portadores. En las prime- 
ras observaciones experimentales, se observaba un comportamiento lineal con dicha 
densidad de portadores, salvo, como hemos mencionado antes, en situación de medio 
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Figura 5.1: Conductividad mínima en el grafeno. En esta figura se muestra el va- 
lor de la conductividad en torno al punto de neutralidad de carga, como función 
de la movilidad de los portadores y, para diferentes muestras. La mayoría de las 
muestras muestran un valor de 0nin en torno al valor 4e?/h, salvo en dos casos (de 
aproximadamente 50). La muestra marcada en color verde mostró inicialmente un 
valor anormalmente alto de Gin, que, tras calentamiento a 400K, adquirió un valor 
próximo al valor medio. Por otra parte, solamente el valor teórico 4e?/1h, predicho 
por la mayor parte de las teorías de transporte en grafeno fue obtenido para una 
muestra. 


llenado n — 0, para la que se alcanza el valor a — 4e?/h. Más tarde, se comprobó 
que el comportamiento lineal de la conductividad con la densidad de portadores no 
es universal en grafeno, dependiendo por ejemplo, de la geometría [55, 56] o del tipo 
de impurezas en la muestra[54]. 

Desde que se publicaron los primeros datos experimentales (e incluso antes) que 
sugerían un valor mínimo universal en ausencia de portadores, se ha realizado un gran 
esfuerzo teórico para tratar de entender dicho fenómeno[57, 58, 56, 59, 60, 61, 62, 63, 
64]. 

Es interesante observar que, en el régimen balístico, es decir, en ausencia de inter- 
acciones y desorden, teóricamente se espera un valor finito y universal de la conduc- 
tividad, en contraposición a lo que se espera para un gas de electrones bidimensional 
usual. El concepto clave para entender este hecho es el llamado “Zitterbewegung” [58]. 


58 Capítulo 5: Desorden topológico y transporte en grafeno. 








o (1/k0) 


B=0 T=10K 











-100 -50 0 50 100 


Figura 5.2: Variación de la conductividad con el potencial de puerta V, en ausencia 
de campo magnético. Se puede apreciar el mínimo en la conductividad Omin a vol- 
tajes de puerta cercanos a cero. El comportamiento en forma de V que presenta la 
figura anterior no es universal, variando de muestra a muestra de un comportamiento 
lineal a sublineal. Este hecho se considera que está relacionado con la naturaleza 
de las impurezas en la muestra[54]. En la figura interior se muestra una imagen de 
microscopía de efecto túnel de una muestra típica de tamaño cercano a la micra. 


Este fenómeno fue propuesto por primera vez por E. Schródinger en los años 30 en 
su análisis de las soluciones en forma de paquetes de onda de la ecuación de Dirac. 
La idea clave subyacente al “Zitterbewegung” es la no conmutatividad del operador 
corriente con el hamiltoniano de Dirac, lo que se traduce en que la corriente no se 
conserva con el tiempo. El mero hecho de la no conservación de la corriente implica 
que las correlaciones entre valores de la corriente a diversos tiempos produzcan un 
valor finito en la conductividad, tal y como se deduce de la fórmula de Kubo. Como 
se puede ver, empleando la fórmula (A.13) del apéndice A, junto con el valor a orden 
árbol del tensor polarización para el grafeno[65), 


E 2 
obtenemos el valor universal Omin = a 


Muy recientemente, han surgido trabajos en los que se considera la posibilidad de 
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que la conductividad mínima no sea universal, dependiendo, entre otros parámetros, 
de la densidad de desorden. Por ejemplo, en la referencia [54] se especula con la 
posibilidad de que la naturaleza de las impurezas sea determinante a la hora de 
determinar el valor de la conductividad mínima y de la forma de la conductividad 
como función del voltaje de puerta. En esta línea de razonamiento, se ha investigado 
con más detalle el papel de las impurezas cargadas situadas en el sustrato, en la 
situación en la que el potencial químico está lejos del punto de Dirac. En este caso, se 
obtiene, tanto teórica como experimentalmente[66, 67] un valor para la conductividad 
proporcional al cociente entre la densidad de portadores n y la densidad de impurezas 
mu rre El escenario anterior se puede generalizar a la situación en la que 
el voltaje externo se sustituye por un sustrato cargado[68] que genera fluctuaciones en 
la densidad de portadores, encontrándose que el valor de la conductividad es de nuevo 
no universal, y verificándose la ley de escala anterior, donde el papel de densidad de 
impurezas lo juega la fluctuación de carga ón inducida por el sustrato. 

Exceptuando el trabajo de Nomura y MacDonald[54], todos los resultados ante- 
riores están basados en la presencia de una densidad de portadores en la muestra, 
ya sea inducida por la aplicación de un voltaje de puerta externo o generada por la 
presencia de impurezas cargadas en el sustrato. Una de las preguntas que nos vamos 
a hacer en este capítulo es qué cabe esperar en cuanto a la universalidad del valor de 
la conductividad mínima en presencia de desorden topológico, donde la densidad de 
portadores es inicialmente nula y no existe un mecanismo externo que induzca den- 
sidad de portadores en el grafeno, tal y como pasa con las impurezas cargadas. Un 
resultado interesante al que llegaremos es que, al menos a nivel semiclásico, el valor 
de la conductividad mínima no es universal, y depende inversamente de la densidad 
de desorden, Y n re comportamiento típico de los sistemas difusivos. 


min ni 


imp: Cmin 


5.2 El modelo. 


En esta sección vamos a recordar brevemente el modelo propuesto para estudiar 
el efecto del desorden topológico sobre la estructura electrónica del grafeno. 

Como hemos visto, un defecto topológico consiste en la substitución de un anillo 
hexagonal por un anillo poligonal de n lados, donde n es distinto de 6. Para definir 
mejor el modelo, nos vamos a centrar en pentágonos y heptágonos, que son los defectos 
que se pueden encontrar con mayor probabilidad en muestras reales. 

En el capítulo anterior vimos que la inclusión de dichos defectos sobre una muestra 
plana de grafeno genera dos tipos de efectos: curvatura local, e introducción de fases 
extra de tipo Aharonov-Bohm. 

Sin despreciar el efecto de dichas fases extra, cuya presencia genera efectos muy 
interesantes sobre la estructura electrónica [9, 30, 33, 69], en el presente capítulo nos 
vamos a centrar en los efectos que la curvatura induce sobre la estructura electrónica 
de la lámina de grafeno curvada, siguiendo el modelo de fermiones de Dirac en espacios 
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curvos descrito en los capítulos anteriores. Vamos a describir la situación realista en 
la que el número de defectos pentagonales y heptagonales es el mismo. Esta condición 
nos permite, primero, asegurar la planaridad asintótica de la muestra, es decir, vamos 
a permitir que la muestra se curve localmente en puntos arbitrarios, pero que a gran 
escala no se separe en manera apreciable de la situación plana, y segundo, hacer que, 
de entrada, una de las fases extra que inducen estos defectos, concretamente la fase 
que da cuenta de la mezcla de los puntos de Fermi o valles asociada a la matriz 7, 
sea nula. 

La segunda fase, relacionada con la no conmutatividad de los operadores de las 
holonomías de traslación y de mezcla de valles, va asociada con la matriz 73 6 1, está 
acompañada por un parámetro y que puede tomar los valores SE y 0. Supondremos 
que dicho parámetro y es nulo. Esto hace que los resultados que obtengamos no sean 
totalmente generales. Sin embargo, la inclusión en el modelo de estos campos gauge 
que no consideramos puede realizarse de forma sencilla, puesto que, al ser formalmente 
similares al campo gauge asociado a la conexión de espín, tras el promedio sobre 
desorden generará un término de interacción de tipo corriente-corriente. 

En el capítulo anterior también vimos la forma de obtener la métrica de un con- 
junto de defectos, modificando la métrica de un único defecto cónico, a través de una 
transformación de coordenadas. En esta ocasión vamos a obtener la métrica a través 
de la dinámica de los defectos en la superficie curva. Esta forma de obtener el tensor 
métrico tiene la ventaja de permitirnos calcular correladores del factor conforme en 
distintos puntos de la superficie. Este tensor métrico, en coordenadas conformes, se 
escribe como: 


Las ecuaciones que describen la dinámica de los defectos se reducen a única ecuación 
para el factor conforme A(r): 


N 


VA(r) => njó(r — ry), (5.2) 


J 
cuya solución en coordenadas conformes es 


od r—Tj; 
A(r) = > > log | (5.3) 
J 








10) 


El parámetro uz está relacionado con el ángulo de defecto o exceso presente en el 
defecto y toma los valores y = Eo, siendo el valor positivo para pentágonos y negativa 
para heptágonos. a* es una constante del orden de la constante de red, que se puede 
interpretar como el radio del núcleo del defecto. El valor de y permite desarrollar 
en serie los términos dependientes de y en el hamiltoniano curvo (4.13) a través de 
(5.3). La dependencia explícita de los diferentes términos que aparecen en (4.13) 
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con el factor conforme A(r) se ha deducido en el capítulo anterior. El hamiltoniano 
resultante es: 


E =iop [14 Ad 0 + FA) (5.4) 


En la ecuación (5.4) aparecen dos términos extra con respecto a la situación plana. El 
primer término representa una modificación local de la velocidad de Fermi inducida 
por la curvatura de la superficie. El segundo representa un campo magnético efectivo 

Estos términos de desorden son dependientes de las posiciones concretas r; de los 
defectos. En el capítulo anterior vimos que la densidad de estados local se puede 
modificar en un sistema con un conjunto pequeño de defectos. En este caso, nos va a 
interesar la situación en la que, en un sistema idealmente infinito, se ha introducido 
un conjunto infinito de defectos, caracterizado por una densidad Nimp. Además, como 
hemos mencionado en la introducción, las propiedades de transporte no deben de- 
pender de la distribución concreta de defectos. En la siguiente sección estudiaremos 
la manera de construir funciones de transporte promediando sobre la distribución de 
defectos. 


5.3 Promedio sobre desorden. 


Definamos la función generatriz %Z de un sistema de campos fermiónicos Y y Y: 
como la integral de camino 


ZIA]= | Die Te S SIGAS, (5.5) 


El valor esperado de un cierto observable Ó puede obtenerse por medio de derivación 
funcional de la función generatriz 5.5: 


n 


O= 57 


l;=01MZ|A, Jl, (5.6) 
donde hemos definido la función generatriz a través de la integral de camino de los 
campos 

donde Hp se definió en (5.4), y en general, O será un funcional de A(r). 

El promedio sobre desorden se puede definir como 


(O)dis = | DATO), (5.7) 


Para darle carácter operacional a la definición (5.7) es necesario tener una acción 
para A(r). En la mayoría de los trabajos sobre desorden, simplemente se considera 
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la función que define el desorden V(r) como una variable estadística con un valor 
promedio nulo (V) ,,, = 0 y una varianza finita, 


(V(OV() q, = Ad(r — 1). 


Toda la información sobre la densidad de desorden está contenida en el parámetro 
adimensional A. En nuestro caso, podemos hacer una interpretación un poco diferente 
del desorden y determinar unas propiedades similares a partir de lo que sabemos de 
la “dinámica” de los defectos topológicos. 

Cuando definimos la función A(r), vimos que satisfacía la ecuación no homogénea 
(5.2), en donde la curvatura S(r) = Pad nj0(r — r;) es el término independiente. La 
solución de (5.2) se puede escribir formalmente mediante la función de Green K(r—r') 


como: 
= Jarstr, rj)K(r —r'). (5.8) 

donde K' es 
K(r —r”) =l0g 





(5.9) 








Esto nos permite definir una acción Ha para A(r) a partir de K(r—r'): 





HE / drd' Ar) K Ur —r)A(r). (5.10) 


esa 
Es importante notar el término adimensional n;mpa* en la definición de Ha, que jugará 
el mismo papel que juega A en la definición estadística de desorden anteriormente 
mencionada. 
En lo que ABU vamos a trabajar en el espacio de momentos, en el cual Hp tiene 
la forma ((dp) = Ep): 


Ho= a (App). (5.11) 


AMimpa? 
Con todo lo anterior podemos postular una varianza para la función A(p) como: 


Mia 
(AD)A(=p)) = == — 
Pp 
Una vez que tenemos definida la operación de promedio sobre desorden, retorne- 
mos a la definición de valor esperado de un observable O como derivada variacional 
de ln Z, y su promedio sobre desorden: 


E De 
E = Hp .12 
(Oi (el .m2) : f Dre aro ozla d. (542) 


La integral funcional es, en general, imposible de hacer, debido a que la función A(r) 
aparece tanto en el numerador como en el denominador. 
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Este problema se puede solventar de tres maneras diferentes: supersimetría|70], 
la técnica de contornos de Keldysh[71] y el truco de las réplicas[72]. 

En el presente capítulo utilizaremos esta última técnica, consistente en modificar 
la función de partición 4% de tal manera que Z|A, J = 0] = 1 pero manteniendo la 
definición (5.12) intacta. 

La base de esta técnica está en la serie de igualdades formales 


oi ERE A 
In Z = lim ple = lim pe 1), 


que permiten escribir (5.12) como 


LEA R 

(O dis = lim 57 (z eS (5.13) 
El resultado formal es que se puede realizar el cálculo de (O) ,,, calculando Z* en 
vez del logaritmo de Z. El cálculo de 2% no conlleva ningún problema, ya que, en 
términos de integral de camino, 


(Je) Jn 


y el promedio sobre desorden se hace ahora en una exponencial cuyo exponente de- 
pende linealmente de A(r) . El nombre de truco de réplicas procede del hecho de 
haber “replicado” la acción R veces en la exponencial de la integral de camino. He- 
mos pasado de trabajar con un conjunto de campos fermiónicos (4, 1) a trabajar con 
R especies de campos (4, 1r),promediados sobre R. Al final del cálculo se tomará 
el límite R = 0. 

El último paso entonces es el cálculo de 


(a ] DAe Ha zR. 


Empleando (5.5) expresada en el espacio de momentos, calculando Z* con (5.4) y 
realizando la integración funcional sobre A(p), obtenemos el siguiente hamiltoniano 
de “interacción”: 


R 
H= Ho +0 aa DANA (na (5.14) 


donde A= 27 1"Nimpa'? y la suma es sobre índices de réplica a y b. Hp es el hamilto- 
niano libre usual 


Ho = J aiebali" — UP yk +i9M)Va, (5.15) 


Tanto en (5.14) como en (5.15) hemos definido dos tipos de campos fermiónicos: 
campos avanzados, para los cuales el parámetro 1 es positivo en (5.15), y campos 
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retardados, donde 1 es negativo. En presencia de estos dos tipos de campos, y para y = 
0 en (5.15), el hamiltoniano (5.14) posee una simetría O(2NW'). Cuando se introduce 
el término asociado a 7, dicha simetría se rompe y pasa a ser O(N) x O(N). Tal y 
como se explica en la referencia [73], esta ruptura explícita de simetría genera modos 
Goldstone, que como veremos más adelante, se interpretan como los modos de difusión 
del sistema. La matriz M se define, por tanto, como 


1.0 
O 610 


En toda la discusión que se sigue, nos centraremos en estudiar lo que sucede al nivel 
de Fermi, es decir, a w=0. 
La función T'(k, k”) tiene la expresión 


El segundo término en T'(k,k'”) procede de la integración funcional del potencial 
vector en (5.4). Dicha integración se ha realizado fácilmente en el espacio de mo- 
mentos, donde O¡A(r) — k¡A(k). El correlador para el campo magnético efectivo 
A; - 0,A([74]) se ha tomado como 


(5.17) 


(As(P)Aj(—P)) = di¡277Mimpa”. 


Así pues, el término correspondiente al campo magnético efectivo genera un 
término de correlación de corto alcance, mientras que el primer término de T'(k, k”), 
procedente del promedio del término que modifica localmente la velocidad de Fermi, 
es un término anisótropo y fuertemente divergente en el límite de dispersión hacia 
delante: k” — k. Esta divergencia procede del correlador (A(p)A(—p)) a momentos 
pequeños, o, equivalentemente, a distancias muy largas. Vamos a regularizar esta 
divergencia en (5.17) por medio de un cutoff ó: 


o a da e 37 :) | (5.18) 





Podemos darle un significado físico a 0. Si en vez de considerar un propagador de 
largo alcance para los defectos que se comporta como Pm consideramos un propagador 
del tipo 
1 
K(p) = ==, 
(p) 248 


que corresponde, en el espacio real, a cambiar el logaritmo en (5.9) por la función de 
Bessel modificada de segunda especie: 


(5.19) 


K(r—r”) = Ko(9 [r — r'|). (5.20) 
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El término dominante en el desarrollo en serie de (5.20) para pequeños valores de d 
es 
K(r — r') = log(6 |r — r'|). (5.21) 


Podemos interpretar entonces 4 como el inverso de la distancia x sobre la cual dos 
defectos topológicos dejan de estar correlacionados. 

Aunque hayamos encontrado un término de interacción de cuatro fermiones tras 
promediar sobre A(p), hay que recordar que estamos tratando con un problema de 
dispersión elástica de electrones, en el cual los módulos de los momentos de los estados 
entrante y saliente tras un proceso de dispersión por una impureza se conservan. En 
virtud de esta consevación de momento, la expresión (5.18) no es función del módulo 
de los momentos, definiendo la diferencia de ángulos de incidencia y salida d = 0, —0y,, 
y Ó?= a encontramos 





E cos?(H/2) E 
19) = A + 62 :) 22) 


Podemos asumir, además, que los momentos involucrados en nuestro problema están 
en un pequeño entorno de los puntos de Dirac, por lo que podemos suponer que 1/k 
es del orden de la longitud de onda A de dichos estados ( o la longitud de correlación, 
si estuvieran localizados). Por lo tanto podemos estimar ó del orden de A /x, siendo 
x comparable con el tamaño del sistema L. Estas hipótesis no son suficientes para 
caracterizar bien ó, siendo por tanto un parámetro incontrolable en nuestra teoría. 
Solamente realizaremos la hipótesis de que x > A. A pesar de ello, veremos que la 
conductividad Drude calculada resulta ser independiente de este parámetro, y por lo 
tanto, bien definida. 
Para concluir esta sección, vamos a descomponer '(H$) en armónicos, 


DO=M Pre" = Y Dacia, (5.23) 


A 


donde, por conveniencia, hemos definido xn(k) = e Esta descomposición es de 
enorme utilidad, ya que nos va a permitir determinar qué canales de la interacción 
(5.14) van a ser los que realmente van a contribuir en el cálculo de la conductividad. 
Además, nos permite comparar nuestro modelo, que tiene múltiples canales de dis- 
persión, con modelos de corto alcance, en los cuales la dispersión ocurre únicamente 
en el canal s. 

Adelantando acontecimientos, vamos a encontrar que solamente los modos con 
n=0 y m= +1 contribuyen al proceso de difusión electrónica. Los valores explícitos 
de los coeficientes TP, para estos valores de n son (P_, =T',,): 


1 1 
lo = 5 5/4+0(0) : Dj =5-2+01(9). (5.24) 
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5.4 Aproximación de campo medio. El modelo y 
no lineal. 


En esta sección vamos a construir la teoría de campos de baja energía que des- 
cribe el comportamiento a escalas muy grandes del sistema, es decir, a muy grandes 
longitudes de onda|75). 

Sustituyendo la descomposición en armónicos en el hamiltoniano de interacción 
(5.14) tenemos: 


R 
E = Y YA) EAN. (5.25) 


con A = 274 4Nimpa*?. Podemos transformar la interacción corriente-corriente (5.25) 
en una interacción de tipo densidad-densidad aplicando una identidad de Fierz[76|: 


(1 Mip2)(13N Ya) = Y (01 B:V4)(V3 MB¿N Yo), (5.26) 


2 


donde M y N son matrices genéricas de SU(2) y B; = ([1, 71,72, ys) es una base de 
matrices de SU (2). 
El hamiltoniano (5.25) se escribe, entonces, como 


R A A a En Em EN 
Hini = 02 Y farra) Y Pr () xk) [(darW)(botba) a (bar51)(Deyspa)] 


(5.27) 
donde hemos definido yg = 24. Es interesante notar que el hamiltoniano de interacción, 
tanto en la forma (5.25) como en (5.27), es invariante bajo transformaciones globales 
continuas U,(1) quirales, en contraste con una interacción de tipo (4,4)? ([77)), 
genuina de un problema de dispersión elástica isótropa de fermiones, que es invariante 
bajo U,(1) discreta. Pese a esta diferencia, veremos que en ambos casos es posible 
encontrar una solución de campo medio que rompa esta simetría, tanto en el caso 
continuo como en el discreto. 

Podemos desacoplar la interacción efectiva de cuatro fermiones(5.27) por medio 
de una transformación de Hubbard-Stratonovitch 


Ha = (at) LT, (03 +1) +1 E Pxa(i0da (Qu — sl) úl. (528) 


Los campos Q% y II% son campos auxiliares matriciales de dimensión 2NV, debido a 
que hemos definido campos fermiónicos retardados y avanzados. Podemos simplificar 
los cálculos que vienen a continuación realizando una transformación unitaria que 
diagonaliza el hamiltoniano libre (5.15) y deja inalterada la medida de integración 
funcional: Y > YU*+, Q, > UQ.U*, IL, > UN, U*, ys > Uy5U*. Sin pérdida de 


generalidad, renombraremos los campos transformados como los originales. 
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Escribiendo la integral de camino o a en función de los campos auxiliares 


Ss _ 
Q*% y TI% e integrando los campos fermiónicos (4, 17), obtenemos la siguiente acción 
efectiva: 





jo JT E TO 2 +18), (5.29) 
donde hemos definido G; como 
G7=G¿'+i9M +10 Taxn(kNQn — ¿nsIa), (5.30) 
y 
a E E ) on (5.31) 


La integral de camino escrita en términos de los campos auxiliares, estará dominada 
por las configuraciones de los campos Q% y 11% cercanos a la solución de campo 
medio, definida a partir de las ecuaciones de punto de silla para dichos campos, 


Hess _ y 9Hes5 
On) (lr) 
Es fácil comprobar que las ecuaciones anteriores pueden agruparse en un único con- 


junto de ecuaciones para los campos compuestos S, = Q, — vy5Hl,,, verificándose que 
52 = S5Sn = Q? + IP?, en virtud de las reglas de conmutación de los campos Q% y 





=0. 


o, SH E 
T1%. La ecuación 7 =0 se escribe como 





S») 
1 Xn(k) 
A O - A 9.32 
2gu2, (Sn) f E +i9M +12 nxn(k) (Sn) E 


Siguiendo la referencia [78], consideraremos como solución de prueba para estas ecua- 
ciones 





donde f es una función a determinar. En el límite y — 0, dicha función tiene la forma 
K 2n -3 
F = pe (0 == 1) ] ] (5.34) 
0 


expresión en la que hemos tenido que introducir un momento de corte ultravioleta 
K. Ahora se pone de manifiesto el hecho de tener una solución que rompa la simetría 
0, (1) continua, ya que la solución (5.33) no fija de manera unívoca las soluciones 
para (Qo) y (Ho). Podemos elegir entonces aquella solución en la que (Qo) = FM y 
(Ho) =0. 

Una vez elegida una solución de campo medio (5.33), el tiempo de vida media o 
tiempo de relajación de una partícula se define como 


a 5.35 
UE (0% 1) , (5.35) 
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resultando ser una constante, y con una dependencia con la densidad de defectos 
típicamente no perturbativa: 77*(g) = exp(1/g). Este valor de la vida media depende 
de ambos regularizadores, K y $. La dependencia con K' puede ser eliminada por 
medio de técnicas de grupo de renormalización [57], técnica que no se puede aplicar 
para eliminar la dependencia con $. Las propiedades de una partícula, como la vida 
media (5.35) y, como veremos a continuación, la densidad de estados, fuertemente 
dependen de dicho parámetro. 

La densidad de estados promedio en el nivel de Fermi puede calcularse a través 
de la expresión 


p(0) = lim qm 1). (5.36) 


G*(k) es la parte retardada ( y = —1) de la función de Green (5.30). En función de 
la vida media, p(0) se escribe como 


211 
- glyut 27" 





p(0) (5.37) 


La densidad de estados finita (5.37) nos va a permitir utilizar la substitución usual 








Vr 























Sy ñ 
0 0 
o o 


Figura 5.3: Densidad de estados e inversa de la vida media como función de la energía. 
El valor de la densidad de estados escala con 91/?, mientras que 1/7 escala con 1/91/?. 


para la integración de momentos en gases de electrones bidimensionales de hamilto- 
niano parabólico: 
2 
ld — plo) f ae , UPkK > Ep. 

Esta substitución tendrá efectos no triviales, ya que, en el caso de tener grafeno 
perfecto, la densidad de estados escala con la energía, como ya sabemos, de la siguiente 
manera polw) — lw|, impidiendo en general utilizar la substitución anterior. 

Antes de calcular las fluctuaciones cuánticas en torno a la solución de simetría 
quiral rota, debemos hacer otro comentario. La ecuación (5.32) admite otra solución, 
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aparte de la definida por las ecuaciones (5.33) y (5.34), que es (S,,) = 0, y que corres- 
ponde a la fase donde la simetría no está rota. La densidad de estados correspondiente 
a esta solución se hace nula en el nivel de Fermi, y eventualmente tendría la forma, 
en función de la energía, de p(w) — lw|”, siendo a: una función de g[60]. 

Para determinar cual de las dos soluciones es de menor energía, calculamos el 
potencial efectivo como función de f. El potencial efectivo V.,f calculado es el mismo 
que el correspondiente al modelo de Nambu-Jona Lasinio [79], tomando la expresa- 


mente (II¿) cero: 
2 je 
epi = laz iS ad - (or (3) - ) (5.38) 


La solución f = 0, a la que AN un valor para el potencial V.;¿ = 0 resulta 


ser un máximo, mientras que la solución f es es efectivamente un mínimo, y co- 
K? 12 : 4 , 
rresponde al valor V.yy = —%7f”. Por lo tanto, tiene sentido calcular las correcciones 
al valor de punto de silla de la integral de camino en torno a la solución (5.34). 
Podemos desarrollar (Q,, y Il, en torno a sus valores de campo medio Q,, = (Q») + 
90, y 5, = (1,,) + 011, en (5.29) y escribir 


Hop A 0 ( Hosp) 7 EI (Quim > Um) E (5.39) 
El símbolo * significa que la derivada variacional de H.¿y está evaluada en los valores 
de los campos en la aproximación de punto de silla. En la ecuación (5.39) aparecerán 
términos cruzados de 9Q,, y ÓlL,,, , esas derivadas funcionales originan términos nulos, 
y por lo tanto los campos anteriores no se acoplan entre ellos. Efectivamente, las 
derivadas de los campos tienen la siguiente forma: 





orion = Er + - ao) (dp) ET Pnxldmo + 90, 
TrG(IG( + 9) m, (5.40) 
y 
sm, ES sy, — aL IT? — N dp)(dq) > TP mxXnl2)xm(p + 9)5H,, 
¿11,011 Am 
Tr ClsCl + q)011 (5.41) 


El término Qi SII, es proporcional a TrG(p)ysG(p + q) que no contribuirá al 
término de correlación de difusión correspondiente (9Q(q)011(—q)). 

Las funciones espectrales se encuentran concentradas en torno a la energía de 
Fermi, por lo que podemos restringir la integración en la variable q a valores en 
torno al nivel de Fermi, q < Kp, y hacer la aproximación xm(p +4) Y Xml(P), 


produciendo un error despreciable mientras la integración sobre k se realice en un 
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entorno muy pequeño alrededor de Kp. En la ecuación (5.40), la integral no nula 
para el producto GFG4, es decir, para el sector no diagonal 9Q*7, da lugar a un 
comportamiento de polo difusivo para los campos 9Q,,. Por el contrario, empleando 
la propiedad de las funciones de Green y¿G*+4y, = —G4*, veremos que, en el caso 
de los modos dII,,, el polo difusivo aparece en el producto G*G* y su compañero 
avanzado-avanzado, indicándonos que el canal difusivo será el canal asociado a las 
partes diagonales JII++»+"" de los campos JII,,. El signo — procendente de la simetría 
de G con ys corrige el signo relativo de los segundos términos en el lado derecho de 
las ecuaciones (5.40) y (5.41). 

Realicemos la integración en k en la expresión (5.40). La misma integración se 
efectuará en (5.41), aplicando la relación entre G* y G%. Definamos la cantidad 
CmnlN, q) como 


Cnnlm, 9) =5 [(dpxta) xp TE") c+ 9) (5.42) 


cone=pn+ >La estructura diagonal de G*+4, nos permite hacer uso de la relación 





A == 0 
| O A Il E ei ñ (5.43) 


1 
—vupPp—ie —upP|p+q|+ie 
Efectuando el producto y la traza, la expresión anterior queda 


1 1 1 1 
TrGRG4 —= E 7 , Ss (5.44) 
UFp—1ieUplp+q|+ie  UFp+icU_ |p + ql — ie 





Inmediatamente vemos que el primer término en el lado derecho de esta igualdad es 
el complejo conjugado del segundo término, y por tanto podemos escribir 





dl 1 
TrGFG* =2Re : =). (5.45) 
vPp — ie ur |p + ql +e 


Como estamos interesados en el comportamiento de largas longitudes de onda de los 
campos 90),,, desarrollaremos (5.45) en torno a pequeños valores de q. 

El primer término de este desarrollo, C%,, consiste en hacer q = 0 en (5.45) y 
hacer la substitución dk —> p(0)dej.. El resultado es 


1 


á 
O e 
mm = P(0) .E 


y (5.46) 


o, desarrollando nuevamente en el límite y — 0: 


O 1 
Cam = 2map" 7 qA(0)óam(27)n + O(1). (5.47) 
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Podemos agrupar el término constante en (5.47) con el término proporcional a 9Q? 
que aparece en (5.29) al desarrollar en torno a la solución de campo medio, originando 
un término de masa para los campos: 


Las q ET » (1) (602 +11). (5.48) 


De este término de masa, vemos que únicamente los modos dQo y dll, son modos no 
masivos ( los modos Goldstone producidos por la ruptura de la simetría O(2N')), y 
que serán los modos colectivos responsables del comportamiento difusivo del sistema, 
de igual manera que ocurre en los gases de electrones bidimensionales[78]. 

El siguiente término en el desarrollo de Cm es 


1 1 i(m+n)0 il 1 
Cam 4 Joy pd (a —idlvpp +10)? * (—vpp — ie) (—upp + 2») 
(5.49) 
Notemos que, en el caso de la dispersión de electrones isótropa y de corto alcance, 
este término es nulo debido a la integración angular. En nuestro caso, sin embargo, la 
presencia de la función e”*"0 permite un término no nulo lineal en q, lo que originará 
un acoplo entre los modos masivos 9Q +1 y el modo sin masa 904 (y lo mismo para 
los modos 011; y dllo). Cambiando la variable de integración y teniendo en cuenta 
que la integración angular es distinta de cero solamente para n= —m+1, CF, puede 
escribirse como 














1 1 
e A o) fa : 5.50 
o, tras hacer y = 0 e integrar por residuos, 
v 
Cm = PO) Ónma Y Sm (27). (5.51) 


El término C?, aunque es el más bajo en el desarrollo en q, solamente da información 
del acoplo entre modos n = +1 y n = 0. Para tener un término dependiente de q 
para ¿Qo, tendremos que calcular el siguiente orden, C?,,. Empleando los métodos 


anteriores, es fácil ver que este término se puede escribir como 





1 1 2 42 
cha = Ine [aan (ato ts 


upp — le | (upp + ie)? M 2upplupp + ie)? 


) (5.52) 


La integración angular de (5.52) nos dice que habrá términos proporcionales a —nm + 2 
junto con términos proporcionales a d,)». Estas condiciones para n y mM generarán 
términos de acoplo entre el modo do y los modos masivos 9(/+2, junto con el término 
buscado para Qp. Al final del cálculo, vamos a integrar los términos masivos en la 
acción (5.29), quedándonos con una acción efectiva para los modos sin masa. La 
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integración de los modos 4Q+> generará un término de orden q*, despreciable frente 
a los términos de orden q?. El resultado para (5.52), teniendo en cuenta únicamente 
el término de dnm es 








nm = Camp (DO (5.53) 


Una vez obtenido el desarrollo para la función Cm podemos escribir la acción total 
para los modos dQ,=0,+1: 


2 2 2 
Hop 7% [o 2 En (E - 1) Qi — Amer T¿9Qo (a + ESO ) Qo 











4 F na Do de 
| DE a (5Q0ÍQ1 + 5Q0éQ +), (5.54) 


e integrar los modos masivos obteniendo, finalmente, 


(QT 
4guz 





JHo= li dasQo (1 + Dg?) 5Qo. (5.55) 


D es la constante de difusión de los modos 4Qp, renormalizada por la interacción 
de dichos modos con los modos masivos 0(/+.La expresión de D en términos de la 
densidad de estados en el nivel de Fermi es 
e. 
212 gp(0)To—T1' 








(5.56) 


A partir de la constante de difusión D) podemos definir el tiempo de transporte como 
Dis DUpT tr. La interpretación física es el tiempo que tarda un electrón inicialmente 
en un estado de momento k en abandonar dicho estado. Esta constante D permitirá 
definir un recorrido libre medio /. Más adelante discutiremos la validez de nuestros 
resultados en función de este recorrido libre medio. 

Finalmente podemos calcular el valor de la conductividad DC en aproximación 
semiclásica por medio de la relación de Einstein (ver apéndices): 


pmp = La. (5.57) 


El factor 4 en (5.57) procede de tener degeneraciones de valle y de espín real. El factor 
2 extra proviene de tener dos canales difusivos, 9Qo y 0IIy. Como factor más relevante, 
la conductividad DC depende inversamente del parámetro de acoplo y = 41 1%a'N imp, 
que contiene información de la densidad de defectos presente en el sistema. 

En conexión con la teoría de los modelos q no lineales, podemos reescalar los 
campos óQ¿ como U = naa $0, y rescribir (5.55) en función de U en el espacio 
real como[80, 57): 


Him =0p0 | CxTrVU(r)VU"(x), (5.58) 
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donde ahora po actúa de parámetro de acoplo de la teoría. Este modelo a eventual- 
mente nos permitiría utilizar toda la maquinaria de grupo de renormalización para 
calcular correcciones cuánticas a dpo[8l, 57]. 

Sin necesidad de entrar el el cálculo de dichas correcciones cuánticas, podemos dis- 
cutir brevemente el resultado obtenido para la conductividad dpc. El resultado más 
llamativo, en el contexto de transporte en grafeno desordenado, es que la conductivi- 
dad, a nivel semiclásico, es inversamente dependiente de la concentración de defectos, 
una cualidad distintiva, sin embargo, de los sistemas difusivos. Un comportamiento 
no universal similar se ha encontrado cuando se ha estudiado el efecto sobre el trans- 
porte del desorden formado por impurezas cargadas distribuidas aleatoriamente en el 
sustrato sobre el que se deposita la muestra[82, 83, 84]. Sin embargo, una diferencia 
crucial entre estos trabajos y el aquí expuesto es que la muestra de grafeno estudiada 
se encuentra o bien fuertemente dopada o bien existe una densidad de portadores no 
nula generada por efecto campo inducido por las impurezas cargadas. En nuestro 
caso, la densidad de estados es generada por el desorden, de igual manera que en las 
referencias [57, 85]. 

Otra característica a destacar del modelo presentado aquí es la fuerte dependencia 
de las propiedades a una partícula, como la vida media 1/7 y la densidad de estados 
p(0), con el parámetro $. La situación aquí es incluso peor que en un gas de electrones 
bidimensional sometido a un campo magnético aleatorio con correlaciones de largo 
alcance discutido en la referencia [78]. Allí, aunque el el tiempo de vida media de las 
cuasipartículas es divergente cuando Í> 0, el tiempo de transporte 7;, es finito. En 
nuestro, caso, Ty también depende de la densidad de estados al nivel de Fermi, pero 
ésta diverge para el límite anterior. Como hemos visto, a pesar de esta divergencia, 
el cálculo semiclásico de la conductividad Drude es independiente de $ y por tanto 
finito, debido a la particular manera con la que escala la constante de difusión D con 
la densidad de estados p(0). 

Hemos definido el parámetro Ó de manera fenomenológica como el cociente entre la 
longitud de correlación A de los estados en el nivel de Fermi, si estuvieran localizados, 
y la distancia característica x. Sin entrar en el cálculo de las correcciones cuánticas 
a la conductividad, 3 es esencialmente no controlable. Podemos, aún así, hacer una 
estimación del rango de aplicabilidad de nuestro resultado. El régimen difusivo está 
caracterizado por el recorrido libre medio, definido como l = vp7;y, que es mayor que 
la longitud de localización A pero más pequeño que el tamaño lineal de la muestra, 
L. Podemos estimar este recorrido libre medio a partir de 7 y rho(0), y asumiendo 
que a* ma y x  L. La estimación para el recorrido libre medio es 





DANS A. 
Ln (4) Ion (5.59) 


imp 


Con esta expresión podemos obtener unas cotas superior e inferior a la densidad de 
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defectos Nimp. Empleando A < [, obtenemos 


4 L 


ES 9125 (5.60) 
De manera similar, la condición l < L ofrece la cota superior: 
4 1 
AS (5.61) 


5.5 Conclusiones, problemas abiertos y trabajo fu- 
turo. 


En el presente capítulo hemos abordado la cuestión de cómo se modifica el trans- 
porte en grafeno en presencia curvatura inducida por defectos topológicos. Basándonos 
en el modelo geométrico descrito en los capítulos anteriores, hemos construido un mo- 
delo sigma no lineal que describe el comportamiento de los modos colectivos en una 
muestra de grafeno en presencia de una densidad de defectos topológicos. 

Debido a la divergencia del correlador para el desorden, el modelo presenta una 
fuerte divergencia infrarroja, que “curamos” introduciendo una distancia característica 
arbitraria que acota la acción de dichos defectos. Debido a esto, todas las propiedades 
a una partícula dependen de esta distancia arbitraria. Pese a esto, el modelo nos ha 
permitido determinar un valor para la conductividad a frecuencia cero que es inde- 
pendiente de esta escala de longitudes arbitraria, haciendo que el modelo sea válido 
cuando se use para estimar cantidades a dos partículas. 

Otro detalle interesante es la dependencia del valor de la conductividad con la 
densidad de desorden. Contrariamente a todos los modelos teóricos no basados en 
impurezas cargadas, el valor encontrado es en esencia no perturbativo, es decir, el 
valor encontrado no es accesible por medio de teoría de perturbaciones en la densidad 
de desorden Mimp, método que nos daría un valor universal. En nuestro caso, el valor 
de la conductividad es el esperado en un sistema difusivo, ya que varía con la inversa 
de la densidad de desorden. En los casos anteriores, la densidad de portadores no nula 
hace que básicamente el grafeno se comporte como un buen conductor, apareciendo 
de forma lógica la física de metales desordenados. En nuestro caso, el sistema ha 
generado una densidad de estados no nula, que si bien depende de la escala arbitraria 
0, es capaz de hacer que, a nivel semiclásico, el sistema se comporte como un sistema 
difusivo. 

El valor de la conductividad a frecuencia cero presentado en este capítulo ha sido 
calculado a nivel semiclásico, lo que hace que no esté justificado compararlo con 
los datos experimentales. La pregunta siguiente que debemos hacernos es cómo las 
fluctuaciones cuánticas afectan a este resultado. 

En trabajos relacionados con el expuesto aquí[63, 64] se apunta a que el modelo 
sigma descrito en la ecuación (5.58) debe ir acompañado de un término topológico 
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( término de Wess-Zumino). La técnica del grupo de renormalización nos permitirá 
estudiar el flujo de la conductividad con la energía, que, según los trabajos antes 
citados, puede converger al valor universal. Esta consideración se basa en el análisis 
de los posibles puntos fijos de la conductividad en función de las simetrías del modelo 
fermiónico con desorden original. Debido al valor concreto del parámetro de acoplo, 
la validez de esta suposición es cuestionable. La investigación de la modificación de 
los puntos fijos en la conductividad se realizará en un futuro. 


Capítulo 6 


Formación de haces polarizados en 
el espín de valle en grafeno. 


6.1 Introducción. 


Como ya se vió en el primer capítulo, la estructura electrónica a baja energía del 
grafeno viene bien descrita por un modelo de enlace apretado a primeros vecinos para 
los orbitales hibridados r en el grafeno. Dicho modelo origina dos bandas de energía 
que se tocan mutuamente en dos puntos no equivalentes ( puntos de Dirac). También 
vimos que estos dos puntos eran el origen de una simetría adicional, que llamamos 
espín de valle o degeneración de valle. En situación de medio llenado, en la que tene- 
mos un electrón por orbital rm, vimos que el nivel de Fermi se encuentra exactamente 
en dichos puntos. También hemos visto que las excitaciones de baja energía dentro 
de esta aproximación tienen estructura espinorial y obedecen la ecuación de Dirac 
sin masa. El origen de esta estructura espinorial en las funciones de onda está en la 
estructura concreta de la red hexagonal que presenta el grafeno y en la degeneración 
de valle. El espín real electrónico daría una estructura SU(2) adicional, sin embargo, 
dicho espín no juega un papel de relevancia en los fenómenos descritos en esta ex- 
posición, por lo que ha sido y sigue siendo considerado aquí como una degeneración 
que únicamente multiplica por dos cualquier resultado que se pueda obtener. Tal y 
como mencionamos en el primer capítulo, esta descripción en términos de espinores 
de Dirac solamente es válida a muy bajas energías, del orden de 0.6 a 0.7 eV, dejando 
de ser lineal por encima de estas energías, según se desprende de medidas de ARPES 
realizadas en muestras de grafito [18]. 

En este capítulo vamos a ver un ejemplo sencillo donde este alejamiento de la des- 
cripción de Dirac tiene efectos observables. Vamos a considerar uno de los dispositivos 
más simples que se pueden diseñar, una unión n — p— n”, donde la segunda región 
n ha sido dopada hasta el punto en el que el nivel de Fermi se sitúe a una energía 
donde la distorsión trigonal es apreciable. Veremos como este simple hecho hace que, 
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Figura 6.1: Superficies de energía constante en torno a los puntos H en grafito. La 
isotropía cesa al alcanzarse energías en torno a 0.6 eV. Figura extraída de la referencia 
[18]. 


bajo unas condiciones concretas, se formen dos haces electrónicos polarizados en el 
grado de libertad de valle. 


Las uniones p — nm son los sistemas más sencillos que se pueden plantear cuando 
se pretende estudiar propiedades de transporte en sistemas cuánticos. En el caso 
del grafeno, resultan ser muy útiles, ya que la transmisión electrónica a través de 
barreras de potencial presenta diferencias sustanciales frente a sistemas descritos por 
el hamiltoniano usual en los gases de electrones bidimensionales. La diferencia más 
relevante es la paradoja de Klein[86). 

Dicho efecto consiste en una transmisión anormalmente alta de electrones a través 
de una barrera de potencial. La transmisión resulta ser uno independientemente de 
la anchura de la barrera, cuando la masa del electrón incidente es despreciable o nula, 
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como ocurre en el grafeno [11]. Sin entrar en el detalle matemático, este efecto es 
fácilmente entendible, si recordamos que los estados electrónicos a baja energía vienen 
descritos por funciones de onda con helicidad bien definida que debe conservarse en el 
proceso, ya que la barrera no conecta estados de distinta helicidad. Concretamente, 
para un estado incidente de momento k% perpendicular a la barrera de potencial, el 
estado reflejado será ortogonal a dicho estado incidente, y el coeficiente de reflexión 
idénticamente nulo. Por lo tanto, la transmisión ha de ser uno, lo que significa que el 
electrón accederá a un estado propagante dentro de la barrera. Esto es exactamente 
lo que ocurre, puesto que dentro de la barrera el electrón encontrará accesibles los 
estados de la banda de valencia, y concretamente, el estado de momento —kí!, el 
cual es inmediato comprobar que dará una transmisión perfecta, de nuevo por mera 
conservación de la helicidad. 

Esta breve discusión sirve para poner de manifiesto que, en el caso del grafeno, 
como ya hemos mencionado antes, es posible encontrar comportamientos en la trans- 
misión de electrones en el grafeno que no son esperables a primera vista en un semi- 
conductor usual. En la siguiente sección vamos a introducir en detalle el concepto de 
distorsión trigonal y veremos cómo afecta a la transmisión de electrones, formando, 
con una propuesta de diseño de unión n—p—n*, un sistema de dos haces electrónicos 
divergentes polarizados en el valle, así como su relación con el concepto recientemente 
acuñado de “valletrónica”. A continuación veremos que ésta no es una propiedad ex- 
clusiva del grafeno, y cómo los cristales fotónicos, de nuevo bajo ciertas condiciones, 
pueden mostrar el mismo comportamiento. 


6.2 Distorsión trigonal, dispositivos y “valletrónica”. 


En esta sección introduciremos el concepto de distorsión trigonal (TW, de sus 
siglas en inglés), mostraremos el cambio cualitativo que la TW produce en la trans- 
misión de electrones por una barrera de potencial, dando algunos números, pensando 
en la posible realización experimental, y acabaremos explicando por qué este efecto 
puede ser relevante, situándolo en el contexto de la “valletrónica”. 

Escribamos, de nuevo, el hamiltoniano en aproximación de enlace apretado para 
los electrones r en el grafeno, incluyendo ambos valles: 


0  óx(k) 0 0 
_| ó(k) 0 y y 
cd 0 0. aa ls 0 
0 0 HK) 0 
siendo p(k) = —t (ass + e73ho cos(Ekya)). En esta ocasión, vamos a desarrollar en 


torno a momentos pequeños (y en torno a ambos puntos de Dirac, y vamos a retener 
los dos primeros términos: 
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0 ko—ik, 0 0 
Ko + ik 0 0 0 
HKk)= uv dd si l 
(k) il 0 0 —kg — ik 

0 0 —Kg + 1ky 0 

0 (kz + ik)? 0 0 

_ Upa | (k, — iky)” 0 0 0 

ñ 0 0 A 
0 0 (lo + ¿k,)? 0 


El primer término es el ya familiar término de Dirac, y el segundo es el término 


de TW. 
la relación de dispersión que se obtiene diagonalizando (6.2) es: 





a? a 
E,(k) = do AL Td sz 3kyk2), (6.3) 


en donde hemos definido el parámetro s = +1 como el índice de valle o espín de valle. 


/ l 





y 








Figura 6.2: Líneas de contorno de la relación de dispersión (6.3) en torno al punto 
de Dirac K. Las líneas de contorno en torno al otro punto de Dirac, K”, se obtienen 
haciendo una reflexión especular sobre el eje k,. 


Como se puede ver en la figura (6.2), el efecto que produce la inclusión de términos 
superiores en el hamiltoniano es el de romper la isotropía en el espacio de momentos 
que presenta el término de Dirac alrededor de cada punto de Dirac, aunque estemos 
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por debajo de la singularidad de Van Hove, y por tanto, la identidad de cada cono se 
mantenga. 

Pasemos ahora a analizar brevemente nuestra unión n—p—nN”. En la figura (6.3), 
puede verse una imagen real y un esquema de un dispositivo de estas características. 
Consiste en una pequeña muestra de grafeno depositado sobre un sustrato, general- 
mente una capa de dióxido de silicio, crecido a partir de silicio dopado, y sobre la 
cual se deposita por litografía un conjunto de electrodos. 


, muestra 
electrodo superior 





grafeno grafeno 





Figura 6.3: Esquema e imagen experimental de un dispositivo electrónico típico de 
grafeno. El electrodo inferior está formado por silicio dopado, sobre el que se ha 
hecho crecer una capa de dióxido de silicio. Sobre el 503 se encuentra la lámina de 
grafeno, y por último, se han depositado el conjunto de electrodos superiorres[87]. 


Los dispositivos como el de la figura (6.3) se basan en el efecto campo. Por medio 
de una diferencia de potencial, llamado voltaje de puerta V,, entre los electrodos 
superiores y el electrodo de silicio dopado, es posible controlar la posición del nivel de 
Fermi en toda la muestra de grafeno. La modificación del nivel de Fermi puede hacerse 
también localmente, permitiendo crear distintas barreras de potencial, tal y como 
mostramos esquemáticamente en la figura (6.4). Por último, aplicando un voltaje de 
polarización V, en los extremos de la muestra, podemos estudiar la transmisión de 
electrones a través de cualquier perfil de potenciales que de diseñe. El valor esperado 
de la corriente que circula por la unión viene expresado en función de la velocidad de 
grupo v, y del coeficiente de transmisión t, por medio de la expresión 


(3) = 5 HH? vo. (6.4) 


El cálculo del coeficiente de transmisión es básicamente un ejercicio de libro de texto: 
conocidas las funciones de onda en cada una de las regiones de la barrera, los coefi- 
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Figura 6.4: Perfil de potenciales en la unión n — p—n” considerada en el texto. 


cientes de transmisión y de reflexión se calculan aplicando la condición de continuidad 
de las funciones de onda en cada frontera de la barrera. Un detalle a destacar es que 
no se requiere una condición de contorno para la derivada de la función de onda, 
puesto que la corriente, como ya sabemos, es un bilineal en la función de ondas, y no 
depende de su derivada, al contrario de lo que ocurre en el caso de tener un hamilto- 
niano cuadrático. En el caso de V, = 0, definido en la figura (6.4), los coeficientes de 
transmisión y de reflexión han sido calculados en [11] y en [88], en aproximación de 
barrera abrupta y suave, respectivamente. Tomar una u otra aproximación, aunque 
modifica cuantitativamente el perfil de los coeficientes de transmisión y de reflexión, 
no modifica sus rasgos característicos cualitativos, como el mencionado al final de la 
introducción del presente capítulo, por lo que adoptaremos aquí la aproximación de 
barrera abrupta, por ser la más sencilla de las dos. 

Sin prestar demasiada atención al coeficiente de transmisión, que no nos va a 
aportar ningún fenómeno nuevo al considerar el término de TW en el hamiltoniano 
(ver apéndices), el comportamiento de la corriente que pasa a través de la barrera 
está controlado por la velocidad de grupo: 


Vas = VkEs(k). (6.5) 


El hecho de que el movimiento del electrón no venga determinado precisamente 
por la dirección del momento adquiere relevancia cuando las condiciones de contorno 
sobre las funciones de onda se traducen a condiciones entre los momentos en las 
diferentes zonas de la barrera. Veamos esto con un poco más de detalle, teniendo 
en cuenta solamente el término de Dirac en todas las zonas de una unión n — p, 
incluyendo a continuación el término de TW. 
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Figura 6.5: (a) Conservación de la componente paralela del momento. Las dos posi- 
bles opciones para el momento transmitido. (b) Cinemática de la transmisión. 


Tal y como se muestra esquemáticamente en la figura (6.5), un electrón en la 
banda de conducción con energía E,, momento k; = |E;| (cos 0;, sin 0;), y velocidad de 
grupo vy = Up (cos 0,, sin 0,) entra en la barrera a un estado de la banda de valencia 
con energía E, = E; — Vo = —vp |k¿]. Por conservación de la componente paralela 
del momento a la barrera, los posibles valores del ángulo de transmisión son 0, y 
0, = T —0,, con cosÓ, = cos 0, pero sin0, = —sin0, (figura 6.5,a). Necesariamente 
la velocidad de grupo del electrón transmitido tiene que apuntar hacia el interior 
de la barrera, lo que fuerza a 0, a pertenercer al segundo cuadrante, y |E,| sin 0, = 
—|Ey|sin0,. Notese que esta derivación es independiente del valle (6.5.a). 

Esta ley de Snell es similar a la que se encuentra en metamateriales zurdos, donde 
el índice de refracción es negativo (aquí, el papel de índices de refracción lo hacen las 
energías: n¿/n; = | Es] / |E,|). Más adelante volveremos sobre esta semejanza con los 
materiales zurdos. 

Cuando se incluye el término de TW, lo primero que hay que tener en cuenta es 
la orientación relativa entre la red hexagonal y la barrera. Ahora que la relación de 
dispersión no es isótropa, el resultado de la conservación de la componente paralela 
de k dará resultados distintos, tal y como se aprecia en la figura (6.6). 

Vamos a considerar la situación en la que |E,| << |E+|. Del capítulo uno, sabemos 
que la primera zona de Brillouin puede tomarse con forma de un hexágono, en cuyos 
vértices se encuentran los puntos de Dirac. Esta zona de Brillouin formará una 
red hexagonal en el espacio de momentos, que se encuentra rotada un ángulo de 


valor 5 con respecto de la red hexagonal en el espacio real. Si nos encontramos en 
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Figura 6.6: Resultado de la conservación del k paralelo cuando el término de TW es 
considerado. (a)y (b): Las dos posibles orientaciones límite entre la red y la barrera. 
(c) y (d): orientaciones entre las superficies de Fermi a un lado y al otro de la barrera 
de pontencial en dichos casos. 


la situación descrita en la figura (6.6.a), se aprecia que para cualquier valor de k; 
permitido, dentro de un mismo valle, la velocidad de grupo correspondiente tendrá 
muy aproximadamente la misma velocidad de grupo vy, y diferirá en dirección con 
respecto a la velocidad de grupo del otro valle, produciéndose de manera efectiva una 
separación de haces polarizados en el índice de valle. 

En el otro caso extremo, (figura(6.6.b)), los momentos permitidos en la zona de la 
barrera presentan un comportamiento muy diferente. Mientras que uno de los valles 
origina un haz totalmente colimado hacia delante, al otro valle le corresponde una 
zona con mayor curvatura, lo que se traduce en un haz mucho más disperso. 

Podemos pensar en la siguiente realización experimental: una muestra de grafeno 
depositada sobre un sustrato de 5102 sobre la que se ha crecido un conjunto de 
electrodos que simulan el perfil de potenciales de la figura (6.4). Los electrones 
son introducidos en el dispositivo a través de una punta de STM, y se transmiten 
por la muestra bajo una diferencia de potencial de polarización, V;. En la figura 
(6.2) se muestran los trazados de rayos correspondientes a las posibles situaciones 
experimentales para ambos casos extremos, cuando las energías son E, = 0.05eV y 
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Figura 6.7: (a): Trazado de rayos correspondiente a las figuras (6.6.a y c) 
. (b): El mismo trazado, correspondiente a la situación de las figuras (6.6.c y d). En 
ambas situaciones se ha elegido |E,| / |E;| = 35. 








E, = 1.75eV. La distancia entre la punta y la barrera de potencial es tal que el foco 
imagen está ligeramente fuera de la barrera. 

Es necesario, con vistas a la posible realización experimental, hacer unas estima- 
ciones numéricas. Como hemos mencionado antes, el término de TW comienza a ser 
relevante a energías del orden de 0.6 — 0.7eV. Sin embargo, las energías típicas al- 
canzadas experimentalmente en medidas de efecto campo son en torno a 0.31eV [89], 
energías a las cuales el término de Dirac predomina. Sin embargo, es posible alcanzar 
valores más altos para el potencial químico modificando algunos parámetros en el 
dispositivo experimental. 

Normalmente, las muestras de grafeno se depositan en sustratos de 5103 de unos 
300nm de espesor. Este valor para el espesor del sustrato no es casual. La manera 
experimental más eficiente hasta la fecha con la que se identifican las muestras de 
grafeno es por contraste óptico[3], es decir, por observación directa a través de un 
microscopio óptico. Para ser identificada, la muestra de grafeno debe tener suficiente 
contraste con el sustrato desnudo. Este contraste, depende esencialmente del espesor 
del sustrato y de la longitud de onda de la radiación electromagnética que se esté 
utilizando, que cae en el rango del visible (entre los 400nm y los 700nm). Esta cota 
en el rango de longitudes de onda hace que el contraste, definido como la variación 
relativa de intensidad luminosa medida en el sustrato con y sin muestra de grafeno, sea 
máximo a espesores de unos 300nm. Sin embargo, teóricamente[3] se ha encontrado 
que existe otro máximo de contraste a un espesor de sustrato de 90 — 100nm, que, 
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junto con un valor típico de voltaje de puerta de V, = 100V, permiten alcanzar 
unos valores para el potencial químico de = 0.6eV, mucho más cerca del valor que 
necesitamos para empezar a notar el término de TW. 


El efecto de separación de haces polarizados en grafeno que hemos descrito antes 
tiene más importancia de la que aparenta. Recordemos una vez más que el grado de 
libertad de valle corresponde a una simetría SU (2) extra, inicialmente presente en to- 
dos los sistemas con fermiones de Dirac no interactuantes en una red ([90, 91]). Hasta 
hace poco, esta simetría jugaba un papel secundario en las propiedades de transporte 
en grafeno, exceptuando quizá el problema de la ausencia de localización débil, donde 
transiciones entre valles debidas a ciertos tipos de desorden([69]) pueden originar una 
ruptura efectiva de la simetría bajo inversión temporal, y la ausencia de magneto- 
resistencia a campos magnéticos pequeños, señal de localización (o antilocalización) 
débil en sistemas electrónicos([92)). 


El espín de valle adquiere protagonismo a partir del trabajo de A. Rycerz, J. 
Tworzydto y C. W.J. Beenakker[1], donde proponen usarlo de manera similar al espín 
real en dispositivos de espintrónica. Esta semejanza pasa por la posibilidad de crear 
y mantener corrientes polarizadas en el espín de valle y la detección eficiente de dicha 
polarización. En este trabajo se propone la creación de un filtro de valle a partir 
de una tira de grafeno donde se ha construido un nanocontacto cuántico, sometido 
al efecto de un voltaje de puerta local. Debido al tamaño finito W de la lámina 
de grafeno, los valores para la componente transversa del momento se encuentran 
cuantizados, con un espaciamiento de niveles A x 1/W, fuera del nanocontacto, y 
9, x 1/W. dentro del nanocontacto, con A. > A. Los estados que pueden atravesar el 

E OE(k) : a . 
nanocontacto son aquellos que tengan una vy = sign (29) definida (según el signo 
del voltaje de polarización aplicado), y el resto es simplemente contar el número de 
estados permitidos con una v¿ dada dentro y fuera del nanocontacto para ambos valles. 
El resultado final es una mayor población de uno de los dos valles respecto del otro, 
y de forma efectiva, la polarización en el espín de valle de la corriente, entendiendo la 
polarización como el cociente entre estados permitidos dentro del nanocontacto con 
una uy dada y estados totales fuera del nanocontacto con dicha uy. La eficiencia de 
esta polarización va decayendo a medida que el nivel de Fermi aumenta, pasando por 
los valores 1,5,5... 

Gracias a este ejemplo, vemos que tener corrientes polarizadas en el espín de valle 
es de relevancia capital si pretendemos usar dispositivos de grafeno basados en la 
“valletrónica”. La desventaja principal de la propuesta de Rycerz y colaboradores 
es experimental: a día de hoy no hay manera ni de obtener tiras de grafeno con 
bordes bien controlados, ni de fabricar nanocontactos como los requeridos en este 
esquema. La propuesta presentada en este capítulo para la generación de haces pola- 
rizados en espín de valle podría suplir dichas dificultades experimentales, al necesitar 
únicamente de una barrera de potencial electrostático, barreras que se pueden obtener 
sin problema hoy día en los laboratorios. 
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6.3 Analogía óptica. 


El fenómeno de haces polarizados en espín no es un fenómeno exclusivo que aparece 
en el grafeno, si lo consideramos como un fenómeno relacionado con la propagación 
de ondas en medios materiales, y nos olvidamos por un momento de la naturaleza de 
dicha onda. Los ingredientes básicos que han dado lugar a la formación de haces son: 
la propagación de una onda a una energía para la cual el medio en el que se transmite 
dicha onda es anisótropo, y la curva de isofrecuencia presenta zonas planas para un 
amplio rango de momentos. Como vamos a ver a continuación, algunos cristales 
fotónicos van a mostrar el mismo fenómeno de formación de haces, en este caso, de 
haces luminosos, exactamente por la misma razón física por la cual los presenta el 
grafeno. 

Un cristal fotónico[93| (PC) es una estructura artificial que posee una constante 
dieléctrica periódica, diseñada para albergar radiación electromagnética de la misma 
forma que un cristal puede albergar y controlar electrones. De la misma manera 
que los electrones en los cristales reales, la propagación de la luz en el seno de un 
cristal fotónico viene gobernada por el teorema de Boch, y, por lo tanto, solamente 
unos vectores de onda estarán permitidos para la luz, dando lugar a la formación de 
bandas en el cristal fotónico[94] (por supuesto, también puede ocurrir que un rango de 
frecuencias no esté permitido para la propagación de ondas, es decir, pueden aparecer 
gaps en la estructura de bandas.) 

De la misma manera que en la propagación de electrones en cristales, la perio- 
dicidad en los cristales fotónicos hacen que los vectores de onda formen un espacio 
finito, la zona de Brillouin. Como ejemplo concreto de analogía óptica con el grafeno 
consideraremos un cristal fotónico bidimensional de estructura de red triangular. En 
la figura (6.8) podemos ver la estructura de bandas para dicho cristal fotónico con 
dos índices de refracción diferentes, n, = 3.4 para las bandas rojas y na = 2.65 para 
las bandas azules. La zona de Brillouin de este cristal fotónico es, al igual que la del 
grafeno, hexagonal, pero con la diferencia de estar ambas zonas de Brillouin rotadas 
una respecto de la otra en el espacio de momentos. Este hecho es debido simplemente 
a la diferente estructura de red que presentan el grafeno y el cristal fotónico que vamos 
a considerar aquí. 

El segundo detalle que nos va a permitir continuar con la analogía óptica es la 
estructura de las bandas en torno a los puntos K y K” de la zona de brillouin. En la 
figura (6.9) se muestra ampliada dicha zona de la banda. De manera similar al grafeno, 
dos bandas de cristal fotónico, a una frecuencia dada, se tocan en un único punto, 
siendo la banda a continuación lineal, es decir, la luz, en el entorno de esos puntos K 
y K' y a unas frecuencias determinadas, se comporta como si fuera libre, pero a una 
velocidad efectiva dada por v, = Vxw(k). Por supuesto, a medida que consideremos 
frecuencias mayores, la aproximación lineal deja de ser válida, apareciendo el efecto 
de distorsión trigonal. 

Ahora que contamos con una estructura de bandas similar, para hacer la analogía 
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Figura 6.8: Estructura de bandas para dos cristales fotónicos de estructura trigonal 
con diferentes valores del índice de refracción. ni = 3.4 para las bandas rojas y 
na = 2.65 para las bandas azules. En la figura se muestra únicamente la estructura 
de bandas para la zona de Brillouin reducida, formada por los puntos M, K' y TI. 
En la zona de Brillouin completa aparecerá otro punto K”, de forma similar a lo que 
ocurre en la zonda de Brillouin del grafeno. 


entre ambos sistemas completa, necesitamos poder simular con cristales fotónicos 
una unión p—nN. A diferencia del caso electrónico, con cristales fotónicos no podemos 
modificar el llenado de las bandas (esto es, el nivel de Fermi) por medio de un voltaje 
de puerta externo, ya que solamente podemos trabajar con frecuencias aisladas para 
la luz. La forma de solucionar este problema es trabajar con una sola frecuencia pero 
considerar dos cristales fotónicos con estructura de bandas similar, cuyos parámetros 
característicos (índice de refracción y parámetro de red) sean diferentes. Éste es 
el motivo por el cual dibujamos en la figura (6.8) dos estructuras de bandas para 
dos cristales fotónicos de parámetros diferentes. En la figura (6.9) se muestra una 
imagen en detalle de la estructura de bandas en torno al punto XK, donde se observa 
la estructura de conos. Aunque en la figura (6.8) no se muestra, tenemos el mismo 
cono en torno al punto K”. 

En las figuras (6.10) y (6.11) se muestra una simulación numérica de elementos 
finitos del módulo del campo eléctrico procedente de una fuente puntual emitiendo luz 
polarizada en el modo transversal eléctrico con longitud de onda A = 4.08apgc en una 
matriz dieléctrica con índice de refracción n, = 3.04 junto con otro cristal fotónico de 
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Figura 6.9: Detalle de la estructura de bandas en torno al punto K para los dos 
cristales fotónicos considerados en el texto. Debido a la diferencia de parámetros entre 
ambos cristales, el punto en el cual se tocan las bandas estará situado a diferentes 
energías en ambos casos, por lo que la velocidad de grupo en ambas situaciones a 
una misma energía será muy diferente, dando lugar a los fenómenos de separación de 
haces o de colimación. 


la misma estructura pero con índice de refracción ligeramente más pequeño na = 2.65 
en ambos casos. 


6.4 Conclusiones. 


En el capítulo 6 hemos estudiado una manera alternativa a la propuesta en la 
referencia [1] para la creación de haces de electrones polarizados en espín. La pro- 
puesta presentada consiste en una unión n —p=—mn” en la que en uno de sus lados 
se ha introducido un dopado por medio de efecto campo tal que el nivel de Fermi 
se encuentra en la zona de energías donde la distorsión trigonal es apreciable. Las 
superficies de Fermi para cada valle, debido a la distorsión trigonal, son diferentes. 
Esto hace que las direcciones de las corrientes transmitidas a través de la barrera 
tengan direcciones distintas para cada valle, produciéndose un efecto de separación 
de haces electrónicos polarizados en el valle. 

El efecto de polarización de valle es muy sensible a la orientación relativa entre la 
red y la barrera de potencial, encontrándose que solamente se da cuando la barrera es 
paralela a la ordenación zig-zag en la red, por lo que la separación de haces electrónicos 
es uno de los dos casos extremos que pueden ocurrir. El otro efecto, asociado a la 
ordenación “armchair”, consiste en la colimación de uno de los dos haces, mientras 
que el otro sufre una dispersión angular. 

Hemos mostrado que las energías en las que la distorsión trigonal comienza a 
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Figura 6.10: Estructura de dos cristales fotónicos que simula una unión p — n. La 
disposición zig-zag de los cilindros produce el efecto de colimación de haces, a dife- 
rencia con el caso electrónico, en donde esta configuración de la red produce el efecto 
de separación de haces. 


ser apreciable están en el límite de las posibilidades experimentales actuales. No 
obstante, variando el espesor del sustrato de Si03, es posible alcanzar valores de 
dopado correspondientes a las energías de la zona de distorsión trigonal. El valor 
del espesor de sustrato es clave a la hora de poder identificar experimentalmente las 
muestras de grafeno por medio de un microscopio óptico, por lo que su variación juega 
en detrimento de la posible realización experimental del dispositivo propuesto, ya 
que, como hemos dicho, para tal valor es complicado alcanzar las energías necesarias, 
mientras que si tal espesor es modificado, aunque sea levemente, se corre el riesgo de 
no poder observar la muestra de grafeno. Sin embargo, existe otro valor para el cual 
el contraste óptico de la muestra de grafeno con el sustrato es óptimo[3], valor que 
permitiría alcanzar dopados suficientes como para ver los efectos aquí descritos sin 
problema. 

Al final del capítulo hemos construido una analogía de la transmisión electrónica 
en grafeno con la transmisión de luz en cristales fotónicos. En tales sistemas es 
posible encontrar puntos de Dirac en las bandas, obteniéndose propiedades de trans- 
porte análogas a las encontradas en grafeno[95, 96, 97]. Entre esas propiedades se 
encuentra la separación de haces electromagnéticos procedentes de distintos valles, y 
la colimación y dispersión angular de haces. La ventaja de la observación de tales 
fenómenos es la ventaja experimental que conlleva, puesto que estas estructuras son 
de uso actual en los laboratorios de óptica. 
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Figura 6.11: Estructura de dos cristales fotónicos que simula una unión p=nN. La dis- 
posición armchair de los cilindros produce el efecto de separación de haces, en el caso 
electrónico, la configuración de este tipo genera él efecto de colimación electrónica. 


Apéndice A 


Teoría de la respuesta lineal. 


A.l El teorema de fluctuación-disipación. 


En general, los experimentos en Física de la Materia Condensada consisten en 
estudiar la respuesta del sistema bajo la acción de perturbaciones externas. En la gran 
mayoría de los casos, si la perturbación es pequeña, la respuesta del sistema será lineal, 
esto es, proporcional a dicha perturbación. Por tanto, el objeto de estudio en esas 
condiciones es el término de proporcionalidad entre perturbación aplicada y respuesta 
medida, que recibe el nombre de susceptibilidad generalizada x. Entre las múltiples 
aplicaciones que posee la teoría de respuesta lineal, una de las más importantes es 
la respuesta del sistema a campos electromagnéticos externos, de donde se extraerá 
información sobre la susceptibilidad de espín o de carga del sistema. 

En este apéndice vamos a introducir de manera muy breve la teoría de la respuesta 
lineal a través de un caso concreto, que sin embargo ha demostrado ser de la mayor 
importancia en el capítulo 3, consistente en la respuesta del sistema a un campo 
eléctrico (o una diferencia de potencial) externo. 

En general, podemos entender la respuesta y (x) del sistema a una perturbación 
externa F(x) como el valor esperado de cierto operador X: 


X (2) = (d(2)Xp(2)). (4.1) 
En esta ocasión, el promedio funcional se va a realizar empleando la acción S = 
So + 0S[F], donde So es la acción del sistema sin perturbar y 


9S|F] = dxFl(aivl(a)X lo), (A.2) 


es el término correspondiente a la perturbación generada por la fuerza generalizada 
F(x). Al final de los cálculos, haremos F' = O para poder identificar el promedio 
funcional en (A.1) con el valor esperado en el vacío de la teoría original, descrita 
únicamente por So: 


X(2) = 5718 2lPllemo (A.3) 
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Vamos a añadir formalmente otro término a la acción total similar en estructura a 
(A.2) asociado a una segunda fuerza generalizada FF”, de tal manera que (A.3) pase a 
ser ahora 


A ') 
X(2) = (0) =-=—— log Z[F, Fl primo. A. 

(2) =(X) = =p y 108 ZlF Pleno (44) 
La diferencia, aunque formal, asegura que el valor esperado (A.4) sea distinto de cero. 
El otro motivo por el cual se ha introducido esta segunda fuerza generalizada es que 
vamos a considerar que X depende débilmente de F”, y es suficiente considerar en 
(A.4) el término lineal en F”: 


Z Ó y | 1 %] 
X(x1) pa =3E(0) (los z17. Fl | Z[F, TT 





2lPllmoP") emo. (45) 


La susceptibilidad generalizada x definida a través de la relación X = yF" tiene la 
expresión siguiente, en términos de derivadas funcionales: 

192 a 
Z[F, FPle=r=0 + ) 


1 0% 
z Sp IF=P=0 


10) == ar 





(A.6) 


Aplicando las derivadas funcionales a 4 y usando (A.4) la expresión de x función del 
operador X' es 


A A A A A 


x(2,1) =(X(2)X (2) - (X(2)(X(2)) = ((X(0) - (X(2))) (X(2) - (X(2")))). 
(A.7) 


La relación entre X(x1) y la fuerza generalizada F(x') puede escibirse entonces como 


(X(0) / dexte, e) Fe). (A.8) 


El el caso de respuesta lineal a un campo electromagnético electromagnético apli- 
cado A, (x) que actuará como fuerza externa F(x), el sistema responderá originando 
una redistribución espaciotemporal de la carga p(1) = (p(1)) y (o) una corriente neta 
j(x) = (¡(u)). Usando una notación generalizada j, = (p,j) podemos escribir (A.8) 
en este caso particular como: 


jule)= [de Kulo,0)A (0). (A9) 
La invariancia gauge introduce una ligadura en la susceptibilidad K,,(x, 2”): 
Aaa. 
Por supuesto, a partir de la definición general (A.6) podemos escribir (4, = 4/) 


$? 
SA, (1)0A,(x') 





Kyle, a) =-27 Z[Alla=o + (4), (A.10) 
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y, por tanto, K,,(x,u') = K,,(x,x'), es decir, K,,(x, x') es simétrico. 

En la teoría de la electrodinámica en medios materiales, la relación (A.9) es ya 
conocida. Si la escribimos en el espacio de momentos, la ecuación integral pasa a ser 
la siguiente ecuación algebraica: 


julw, q) + Kylw, q)Ay(w, q). (A.11) 


Es natural considerar la situación en la que las fluctuaciones medidas sean de tamaño 
macroscópico, comparable al tamaño macroscópico. Además, experimentalmente es 
usual emplear campos eléctricos y magnéticos cuyas variaciones espaciales sean muy 
pequeñas comparadas con el tamaño típico de las muestras. Esto hace que sea lógico 
tomar el límite q > 0 en (A.11). Empleando la condición de gauge Vó = 0 junto con 
E(w) =wA y podemos reescribir la parte espacial de j,, como 


j(w) = 0(1)E(w), (A.12) 


que es la ley de Ohm en su forma vectorial. De esta forma, queda definido el tensor 
conductividad eléctrica a partir de la susceptibilidad K.,,, 


o(w)= lim lora (A.13) 


q4>0 yw 


A.2 La relación de Einstein. 


La relación de Einstein nos va a permitir calcular la conductividad a frecuencia o 
voltaje cero en una aproximación semiclásica a partir de la densidad de estados del 
sistema en cuestión y la constante de difusión sin necesidad de aplicar la expresión 
(A.13) directamente. 

Si consideramos en general un gas de electrones sometido a un potencial eléctrico 
p(r) a temperatura constante, una variación local de dicho potencial eléctrico origi- 
nará una variación en la densidad de electrones n según la expresión ( en equilibrio, 
el potencial químico cumple y(r) = ep(r)) 


ón(r) = Zelr=es edp(r). (A.14) 


Por otra parte, una variación espacial de la densidad de electrones produce un movi- 
miento en forma de corriente de difusión, que, en situación de equilibrio, tratará de 
contrarrestar dicha variación dP(r): 


j=-—eDVn(r) = De olas VO) (A.15) 
Igualando esta expresión para la corriente con (A.12) y empleando E(r) = —Vo(r) 
llegamos inmediatamente a 
d 
g=2t (A.16) 


du lacas 
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Es bien conocido que la densidad de electrones se puede obtener a partir de la densidad 
de estados sabiendo que 


n= , / “o(ENAE. (A.17) 


Simplemente derivando esta expresión con respecto a y y particularizando en el nivel 
de Fermi, obtenemos la versión final de la relación de Einstein: 
E 


h 


e] p(Ep)D. (4.18) 
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